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ABSTRACT 

The first direct detection of gravitational waves (GW) by the ground-based Advanced LIGO/Virgo 
interferometers is expected to occur within the next few years. These interferometers are designed 
to detect the mergers of compact object binaries composed of neutron stars and/or black holes to 
a fiducial distance of ~ 200 Mpc and a localization region of ~ 100 deg 2 . To maximize the science 
gains from such GW detections it is essential to identify electromagnetic (EM) counterparts. Among 
the wide range of proposed counterparts, the most promising is optical/IR emission powered by the 
radioactive decay of r -process elements synthesized in the neutron-rich merger ejecta - a “kilonova”. 

Here we present detailed simulated observations that encompass a range of strategies for kilonova 
searches during GW follow-up. We utilize these simulations to assess both the detectability of kilonovae 
and our ability to distinguish them from a wide range of contaminating transients in the large GW 
localization regions. We find that if pre-existing deep template images for the GW localization region 
are available, then nightly observations to a depth of i « 24 mag and z ^ 23 mag are required to 
achieve a 95% detection rate; observations that commence within ~ 12 hours of trigger will also 
capture the kilonova peak and provide stronger constraints on the ejecta properties. We also find that 
kilonovae can be robustly separated from other known and hypothetical types of transients utilizing 
cuts on color (i — z > 0.3 mag) and rise time (t T i se < 4 days). In the absence of a pre-existing 
template the observations must reach ~ 1 mag deeper to achieve the same kilonova detection rate, 
but robust rejection of contaminants can still be achieved. Motivated by the results of our simulations 
we discuss the expected performance of current and future wide-held telescopes in achieving these 
observational goals, and find that prior to LSST the Dark Energy Camera (DECam) on the Blanco 
4-m telescope and Hyper Suprime-Cam (HSC) on the Subaru 8-m telescope offer the best kilonova 
discovery potential. 

Subject headings: gamma rays: bursts — gravitational waves: binaries, follow-up searches — methods: 
wide-held searches 


1. INTRODUCTION 

The era of gravitational wave (GW) astronomy is fast 
approaching with the Advanced Laser Interferometer 
Gravitational Wave Observatory (ALIGO, Abbott et al. 
2009) and advanced Virgo (AVirgo, Acernese et al. 2009) 
interferometers expected to make the hrst direct GW de¬ 
tections in the 2015 — 18 time frame. These detectors 
will observe the GW signature resulting from the inspi¬ 
ral and merger of compact objects (e.g., neutron star- 
neutron star (NS-NS), solar mass black hole-black-hole 
(BH-BH), or NS-BH binaries). At design sensitivity, the 
ALIGO/AVIRGO network will be able to detect NS-NS 
mergers to an average distance of ~ 200 Mpc (Abadie et 
al. 2010) and with a typical localization region of ~ 100 
deg 2 (Fairhurst 2009, Nissanke et al. 2011, Aasi et al. 
2014). While this will be a remarkable achievement for 
both physics and astrophysics, the detection of an elec¬ 
tromagnetic (EM) counterpart is essential to maximize 
the science obtained from a GW event. The identifica¬ 
tion of an EM counterpart will provide numerous benehts 
including: improved localization leading to host galaxy 
identihcation, determination of the distance and energy 
scales, identification of the progenitor local environment, 
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and insight into the hydrodynamics of the merger (e.g., 
Sylvestre 2003, Stubbs 2008, Phinney 2009, Stamatikos 
et al. 2009, Fong et al. 2010, Metzger & Berger 2012). 
Furthermore, the detection of an EM counterpart will as¬ 
sist in breaking modeling degeneracies to obtain more ac¬ 
curate measurements of binary properties from the GW 
signal (Hughes & Holz 2003). 

Extensive theoretical and observational work has been 
performed to understand the potential EM emission from 
a merging compact object binary. The range of poten¬ 
tial counterparts covers the entire EM spectrum. The 
most popular counterpart is a short-hard gamma-ray 
burst (SGRB; e.g., Paczynski 1986, Narayan et al. 1992, 
Berger 2014) powered by rapid accretion onto the com¬ 
pact remnant which leads to a strongly beamed rel¬ 
ativistic jet. The evidence for a connection between 
SGRBs and compact object mergers includes the associ¬ 
ation with older stellar populations than for long GRBs 
(Berger et al. 2005, Bloom et al. 2006, Berger 2011, Fong 
et al. 2011, Fong et al. 2013), large physical offsets be¬ 
tween the burst locations and their host galaxies indica¬ 
tive of progenitor kicks (Fong et al. 2010, Fong & Berger 
2013), weak association with the underlying stellar light 
of their hosts (Fong et al. 2010, Fong & Berger 2013), 
and a lack of association with supernovae (Fox et al. 
2005, Soderberg et al. 2006, Berger, Fong & Chornock 
2013, Berger 2014). The relativistic jet can also be ob- 
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served in the x-ray, optical, and radio as non-thermal 
afterglow emission arises from interaction with the ambi¬ 
ent medium. This emission is also affected by relativistic 
beaming but may still be seen off-axis as the jet deceler¬ 
ates and expands. In particular, for large off-axis angles 
(6 > 20j) the emission will be most pronounced in the 
radio on timescales of weeks to years (Nakar & Piran 
2011, Metzger & Berger 2012). 

The merger of a compact object binary involving at 
least one NS is also expected to produce mildly relativis¬ 
tic neutron-rich ejecta, with a typical mass of M e j ~ 
10 -3 — 10 -1 M 0 and a velocity of /3 e j ^ 0.1 — 0.3 (e.g., 
Rosswog et al. 1999, Rosswog 2005, Bauswein et al. 
2013a). At later times, less neutron-rich outflows may 
arise from the resulting accretion disk, but the mass 
and composition will depend sensitively on the nature 
of the remnant (e.g., Metzger et al. 2009a, Dessart 
et al. 2009, Fernandez & Metzger 2013, Fernandez et 
al. 2014). The dynamical ejecta is of particular inter¬ 
est when studying potential EM counterparts since it 
is expected to produce an isotropic quasi-thermal tran¬ 
sient powered by the radioactive decay of heavy elements 
produced via r-process nucleosynthesis in the neutron- 
rich merger ejecta (Li & Paczynski 1998, Rosswog 2005, 
Metzger et al. 2010, Grossman et al. 2014, Tanaka 
et al. 2014). The low ejecta mass and high ejecta ve¬ 
locity combined with the resulting r-process opacities 
(ft ~ 10 — 100 cm 2 g -1 for A ~ 0.3 — 3 /im, Kasen, Bad- 
nell, & Barnes 2013) produces a significant suppression 
of bluer optical emission, resulting in a faint (L ~ 10 40 
erg s -1 at peak), rapid (t ~ a few days), and extremely 
red transient with a peak at ^ 1.5 microns (e.g., Kasen, 
Badnell, & Barnes 2013, Barnes & Kasen 2013, Tanaka 
& Hotokezaka 2013), dubbed a “kilonova.” On longer 
timescales (e.g., ^ years), the ejecta will decelerate via 
interaction with the ambient medium producing long- 
lasting isotropic emission in the radio (Nakar & Piran 
2011 ). 

To investigate the detectability of kilonovae, and 
to asses potential contamination in deep high-cadence, 
wide-held searches for such sources, we carried out a set 
of observations using the Dark Energy Camera (DECam) 
mounted on the CTIO/Blanco 4m telescope (PI: Berger). 
This GW follow-up “dry-run” covered an area typical for 
a GW error region (~ 70 deg 2 ), with two visits per night 
for five nights to a depth of « 24 AB mag in both i- and 
z-band. This rich data set will serve as the framework 
to study the nature of foreground/background contami¬ 
nation and to test the effectiveness of various observing 
strategies and contaminant rejection strategies when per¬ 
forming follow-up observations of a GW event. 

To help provide a framework for these observations, 
and to more broadly investigate the coupled issues of 
kilonova detectability and contaminant rejection, we 
present here a systematic study of the feasibility of kilo¬ 
nova detections and assessing the effects of a population 
of contaminant sources during follow-up of a GW trigger 
using wide-held optical telescopes. This study will guide 
the development of observational strategies for wide-held 
searches that maximize the likelihood of detecting a GW 
optical counterpart. The paper is organized as follows: 
In Section [2] we describe the kilonova models utilized in 
this work and develop the list of contaminating tran¬ 


sients considered in this work. In Section 0 we describe 
observational features that may be used to isolate kilo- 
novae from the population of potential contaminants. In 
Section [4] we present detailed Monte Carlo simulations 
of observations to investigate the detectability of kilono¬ 
vae and contaminant rejection for various observational 
strategies. In Section [5] we address issues that arise 
from the possible lack of a deep pre-existing template. 
We then discuss the effects of modihcations to the stan¬ 
dard kilonovae models on our results in Section 1 In 
Section [7j we asses the expected performance of current 
and future wide-held telescopes in the framework of kilo¬ 
nova detectability and contaminant rejection and develop 
strategies for GW follow-up. Conclusions are presented 
in Section 0 

Throughout this work all magnitudes are given in the 
AB system unless otherwise specihed. For cosmological 
calculations, the Hubble constant is taken to be Hq = 
67.77 km s -1 Mpc -1 (Planck Collaboration et al. 2014). 

2. KILONOVAE AND CONTAMINANTS 

To assess contamination from unrelated transients dur¬ 
ing a wide-held search for a GW optical counterpart, 
we consider various fast transients that may present 
kilonova-like behavior, primarily rapid evolution in their 
light curve, in rise and decline. We also consider Type 
la supernovae, as they are the most common SN type in 
large photometric surveys. 

2.1. Kilonovae 

As discussed in Section jl] the merger of NS-NS or 
NS-BH binaries is expectedto produce a faint, rapidly- 
evolving optical/NIR transient powered by the radioac¬ 
tive decay of r-process elements formed in the ejecta. 
From a wide range of simulations it has been shown 
that the ejecta have a typical mass range of M e j ^ 
10 -3 — 10 -1 Mq with velocities of /3 e j ~ 0.1 — 0.3, 
which depend on the initial mass asymmetry in the bi¬ 
nary (Rosswog et al. 1999, Rosswog 2005, Bauswein et 
al. 2013a, Rosswog, Piran, & Nakar 2013). The ejecta 
are expected to be neutron-rich, with typical electron 
fractions of Y e < 0.1 (Metzger et al. 2010). 

A critical property for understanding kilonova emission 
is the opacity, as the ability of photons to diffuse through 
the ejecta determines the timescale, peak luminosity, and 
the spectral energy distribution of the transient. Met¬ 
zger et al. (2010), using the models of Li & Paczynski 
(1998, LP98 hereafter), found that a typical Fe-peak el¬ 
ement opacity (ft ^ 0.1 cm 2 g -1 for A ~ 0.3 — 1 /im) 
leads to a peak time of t p ~ 1 day, with a bolomet- 
ric peak luminosity of L p ~ 10 41 erg s _1 , and a pho- 
tospheric temperature of T p « 10 4 K (i.e., a peak in 
the UV). However, subsequent work by Kasen, Bad¬ 
nell, & Barnes (2013) and Barnes & Kasen (2013) found 
that r-process material have much larger opacities with 
ft ~ 10 — 100 cm 2 g -1 for A ^ 0.3 — 3 /im. This leads 
to a transient that is longer lived ( t p ~ days), fainter 
(L p ~ 10 40 erg s -1 ), and much redder (T p « 3000 K, 
e.g., a peak in the NIR). These results have been con¬ 
firmed in simulations by other groups (e.g., Tanaka & 
Hotokezaka 2013, Grossman et al. 2014, Tanaka et al. 
2014). 

In this work, we construct kilonova light curves 
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1.2 M 0 O-Ne WD with a 3 M 0 BH. The models of 
Metzger (2012) only include bolometric light curves. To 
produce i- and z-band light curves we compute an ef¬ 
fective temperature using the LP98 models with M e j « 
0.5 M 0 , and v e j ~ 2.8 x 10 4 km s _1 for scenario 1, and 
M e j ~ 1 Mq, and v e - } ~ 4.6 x 10 4 km s -1 for scenario 
2 (see Figure 7 of Metzger 2012). We integrate the re¬ 
sulting blackbody curve over the DECam filters using 
Equation [l] and compare to the bolometric luminosity to 
generate bolometric corrections. For the WD-NS merger, 
we find T p ~ 6000 K leading to mboi — ~ —1.9 mag 

and mboi — m z ~ —2.2 mag. The WD-BH merger calcu¬ 
lation yields T p « 7850 K leading to mboi — ~ —2.2 

mag and mboi — m z ~ —2.6 mag. We assume that these 
corrections are constant in time. The bolometric correc¬ 
tions are applied to the light curves of Metzger (2012) to 
generate the light curves used in this analysis (see Fig- 
ure|2]). 


Fig. 1.— Top: Kilonova spectra from Barnes & Kasen (2013) at 
~ 1 d and ~ 1 week with ejecta properties of M e j = 10 -2 Mq and 
j3 e j = 0.2. Bottom: Light curve for the same model parameters 
computed at a fiducial distance of 200 Mpc in the DECam r, i- 
and z-band filters. 


using the time-resolved spectra of Barnes & Kasen 
(2013). They consider a grid of models with M e j = 
10 -3 ,10 -2 , and 10 -1 M 0 and with /3 e j = 0.1, 0.2, and 
0.3. Unless otherwise noted we only consider the average 
ejecta properties of M e j = 10 -2 M© and /3 e j = 0.2 since 
we do not know the distribution of ejecta properties a 
priori. We comput e lig ht curves using the DECam filter 


transmission curver 
filter 0, at a time t, 


with the luminosity through a 
eing given by 




f d\L\(t, A)0(A) 
f d\(/)(X) 


(i) 


where L\(t , A) is the model specific luminosity and 0(A) 
is the filter transmission as a function of wavelength. A 
representative plot of the spectra at 1 day and 1 week, 
along with the resulting r-, z-, and z-band light curves 
are shown in Figure []] 


2.2. WD-NS/WD-BH Mergers 

The tidal disruption of a WD by a NS or BH bi¬ 
nary companion can produce an isotropic optical tran¬ 
sient powered by the radioactive decay of 56 Ni synthe¬ 
sized in the ejecta. The ejecta have an expected mass of 
M e j ^ 0.3 — 1.0 Mq, velocity v e - } « (1 — 5) x 10 4 km s -1 , 
and 56 Ni mass of Mse Ni « 10 -4,5 — 10 -2,5 M© (Met¬ 
zger 2012). These properties are primarily determined 
by the mass of the WD, with more massive progeni¬ 
tors producing higher ejecta mass, velocities, and 56 Ni 
masses. The resulting light curves have a typical peak 
timescale of t p ~ week and a peak bolometric luminosity 
of L p ~ 10 39 — 10 41-5 erg s -1 . 

The two scenarios we consider here are the merger of: 
(1) a 0.6 Mq C-0 WD with a 1.4 M e NS, and (2) a 


3 http://www.ctio.noao.edu/noao/sites/default/files/ 
DECam/DECam_filters_transmission.txt 

We note that the DECam transmission filters are fairly stan¬ 
dard so these results are relevant for other instruments as well. 


2.3. White Dwarf Accretion Induced Collapse 

The accretion of material onto an O-Ne WD near the 
Chandrasekhar mass can cause the degenerate star to 
undergo electron capture in its interior resulting in the 
formation of a NS, rather than a thermonuclear explo¬ 
sion. This process is known as accretion induced collapse 
(AIC). If the WD is rapidly rotating then a stable ac¬ 
cretion disk will form. The disk will cool and become 
radiatively inefficient, driving powerful outflows with a 
characteristic ejecta mass of M e j ^ 10 -2 M 0 and ve¬ 
locity /3 ej 0.1. The outflowing material has an elec¬ 
tron fraction of Y e ~ 0.5 due to neutrino radiation from 
the proto-NS and synthesizes a few x 10 -3 M 0 of 56 Ni 
powering an optical transient (Metzger et al. 2009b). 
Radiative transfer calculations by Darbha et al. (2010), 
motivated by the models of Metzger et al. (2009b), pro¬ 
vided light curves of AIC events. They found that the 
resulting transient is rapidly evolving with the light curve 
peaking at t p ~ 1 day with a bolometric luminosity of 
L p ~ 2 x 10 41 erg s -1 . 

We use the models of Dabhra et al. (2010) to compute 
i and z-band light curves using the same methodology 
as in Section 2.2 The effective blackbody temperature 


has an assumed fiducial value of T p « 5800 K from the 
models of Dabhra et al. (2010, see their Figure 3). The 
resulting bolometric corrections are therefore m\> 0 \—m\ ^ 
— 1.9 mag and rabol — rn z « —2.2 mag. We assume that 
these corrections are constant in time and apply them 
to the bolometric light curves to generate i- and z-band 
light curves (Figure 2). 


2.4. White Dwarf Double Detonation 

Detonation of a sufficiently large accreted He envelope 
(Mn e ^0.2 Mq) on the surface of a C-0 WD can gener¬ 
ate shock waves which then propagate to the core lead¬ 
ing to the total detonation of the WD. In simulations 
performed by Sim et al. (2012), two possible detona¬ 
tion scenarios were considered. The first, edge-lit double 
detonation (ELDD), occurs when the initial detonation 
takes place at the surface of the WD. This produces a 
transient that is expected to reach a peak bolometric 
luminosity of L p ~ 10 42 erg s -1 at t p ~ 8 d. The sec¬ 
ond scenario is the compression-shock double detonation 
(CSDD), where the core is compressed by the shock wave 
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and C-0 detonation occurs closer to the center of the 
WD. The transient is brighter than in the ELDD sce¬ 
nario with a peak bolometric luminosity of L p ~ 4 x 10 42 
erg s _1 , and a correspondingly longer rise time (t p > 10 
d). In the CSDD scenario the transient is also longer- 
lived with no significant decline in brightness for ^ 30 
d. 

The slow evolution of the CSDD model light curve rel¬ 
ative to kilonovae means that it can be easily removed 
from consideration as a contaminant on the basis of 
timescales alone. However, we still consider the ELDD 
scenario as the more rapid rise time and comparably 
faster decline means that it could serve as a potential con¬ 
taminant. We make use of bolometric light curves from 
the simulations of Sim et al. (2012) to c omp ute DECam 
i- and z-band light curves as in Section [212) The ELDD 
spectrum is approximated as a blackbody with a char¬ 
acteristic temperature of T p « 4800 K using the model 
spectra of Sim et al. (2012, see their Figure 7). The 
bolometric corrections are found to be mboi — un\ « —1.7 
mag and 777+01 — rn z « —1.9 mag. We assume that these 
corrections are constant in time and apply them to the 
bolometric light curves to generate i- and 2 +band light 
curves (Figure 2). 

2.5. Type .la Supernovae 

Another possible WD-based transient is the detonation 
of an outer shell of helium on the surface of the WD. 
This can occur when the mass accretion rate from the 
donor is low enough that burning in the He shell becomes 
thermally unstable. The specific term “type .la” was 
coined as a reference to He detonation in an AM Canum 
Venaticorum system (AM CVns) where the donor star 
is degenerate or semi-degenerate and accreting He rich 
material onto a C-0 or O-Ne WD (Bildsten et al. 2007, 
Shen et al. 2010). The typical envelope mass is M env ^ 
few x 10 -2 — 10“ 1 M 0 and the detonation ejects the 
bulk of the envelope at a typical velocity of v e - } ~ 10 4 
km s -1 . This leads to a transient with a typical peak 
time of t p ~ 2 — 10 d and a peak bolometric luminosity 
of L p ~ 5 x 10 41 — 5 x 10 42 erg s -1 (Shen et al. 2010). 

We produce type .la light curves from the time-resolved 
spectra computed from the models of Shen et al (2010). 
We consider the models 0.6 ± 0.2 M 0 , 0.6 ± 0.3 M 0 , 
1.0 ± 0.05 M 0 , and 1.2 ± 0.05 M 0 , where the first num¬ 
ber is the WD mass and the second number is the mass of 
the He envelope. We integrate the spectra over the DE¬ 
Cam filter transmission curves as described in Section l2Tl 
(Figure [2]). 

2.6. Pan-STARRS Fast Transients 

Drout et al. (2014) discovered ten rapidly evolving 
transients in the Pan-STARRS 1 Medium Deep Survey. 
These events are characterized by a time above half¬ 
maximum of less than 12 days and have typical abso¬ 
lute magnitudes of M « —16.5 to — 20 mag. These ob¬ 
jects have blue colors with g — r < -0.2 at peak (Drout 
et al. 2014). While the existing sample is small these 
events may end up being a significant contaminant in the 
fast cadence, wide-held GW searches investigated here. 
Therefore, we consider the i- and z- band light curves of 
three objects from Drout et al. (2014): PSl-lOah, PS1- 
lObjp, and PSl-llqr. They possess sufficient pre-peak 
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Fig. 2. — Top: i-band light curves of the contaminants considered 
in this work at a fidu cial d i stan ce of 200 Mpc. Bottom: z-band light 
curves. See Sections E2J-E21 for details. 


data for our analysis while providing a representative 
range of timescales and luminosities for the sample (Fig¬ 
ure |2|. 


2.7. Type la Supernovae 


Type la supernova result from the thermonuclear ex¬ 
plosion of a C-0 WD. While these SNe have timescales 
much longer than those of kilonovae and other transients 
considered in this section, their high volumetric rate and 
large luminosity mean that a deep optical search for kilo¬ 
novae will detect large numbers of Type la SNe. This 
makes them a particularly likely contaminant in such 
searches. We focus our analysis on the event SN2011fe 
discovered in M101 by the Palomar Transit Factory (Nu¬ 
gent et al. 2011a) due to the high quality of available 
data. In addition to being one of the closest known 
Type la SNe (/i = 29.04 ±0.19, Shappee & Stanek 2011), 
SN2011fe was discovered a mere 12 hours after explosion 
(Nugent et al. 2011b) leading to coverage starting at -15 
days relative to the R-band peak (Pereira et al. 2013). 
In this work, we use the 32 spectral observations com¬ 
piled by Pereira et al. (2013) to produce DECam i- and 
z-band light curves as a function of redshift, inclu ding 


K-corrections, using the process described in Section |2T 
(Figure [2| 
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TABLE 1 

Expected Rates For Fast Transients 


Object 

^vol 

2max 

7^-area 

A 

Reference 


(Mpc -3 yr -1 ) 


(deg -2 yr -1 ) 



Type .Ia 

10“ 6 

0.5 

0.7 

1.4 

Bildsten et al. 2007 

WD-NS mergers 

10“ 5 

0.06 

2 x 10 -2 

3 x 10 -2 

Thompson 2009 

WD-BH mergers 

10“ 5 

0.06 

2 x 10 -2 

3 x 10 -2 

Fryer et al.1999 

AIC 

10“ 6 

0.06 

2 x 10 -3 

3 x 10 -3 

Darbha et al. 2010 

ELDD 

3xl0 -7 

0.14 

6 x 10 -3 

1 x 10 -2 


Pan-STARRS fast 

5xl0 -6 

0.5 

3.5 

7 

Drout et al. 2014 


Note. — Expected rates for the various contaminants considered in Section [ 2 ] 1Z area is computed assuming an isotropic distribution 
of sources in a volume defined by t he co moving volume at z max . The column Arefers to the number of events expected during a search 
covering 100 deg 2 for 7 days. See § |2.10| or details. 


2.8. Stellar Flares 

Stellar flares, primarily from M dwarfs, are another 
potential contaminant in fast cadence wide-held optical 
searches (Becker et al. 2004, Kulkarni & Rau 2006, 
Berger et al. 2013). These hares have characteristic 
timescales of minutes to hours (e.g., Berger et al. 2013) 
and temperature of T p ~ 10 4 K, leading to an appearance 
as fast and blue transients. As a result, M dwarf hares 
will only appear as a signihcant contaminant in kilonova 
searches that rely on single-epoch detection, primarily in 
g- and r-band in which the kilonova timescale is shorter 
than in the redder bands. We return to this point in 
Section [ 7 ] The sky-projected M dwarf hare rate at the 
depths relevant for kilonova searches is large, ^ 10 deg -2 
day -1 (Kulkarni & Rau 2006). On the other hand, as was 
demonstrated by Berger et al. (2013), deep observations 
in the i- and z-band can easily reveal the underlying M 
dwarfs and hence reduce the potential for contamination 
from hares. 


2.9. AGN Variability 

We also consider short-duration optical variability 
from active galactic nuclei (AGNs). AGN typically show 
strong variability on month or year-long timescales, mak¬ 
ing AGN outbursts too long in duration to be confused 
with a kilonovae. However, recent work has detected 
rare cases of day-long or shorter variability in the NIR. 
Genzel et al. (2003) observed an H-band hare in Sgr A* 
during which the source brightened by a factor of 5. This 
hare lasted ^ 30 min, with a characteristic rise time of 
^2 — 5 min. Such short timescale hares can be ruled 
out as kilonova contaminants by requiring a detection on 
more than a single night. Flares with longer durations 
are also unlikely to be signihcant contaminants. Totani 
et al. (2005) performed V- and z-band observations of six 
AGN, selected on the basis of their variability, over a pe¬ 
riod of four days. They found typical variability of only 
^ 1 — 5% above their host galaxy luminosity. More recent 
work by Ruan et al. (2012) studied the variability of a 
sample of Fermi selected blazars. They found that these 
sources showed a typical rms variability of ~ 0.6 mag on 
timescales of ^ 10 d. Therefore, AGN are unlikely to 
exhibit sufficiently large variability to be confused with 
a kilonova. 


2.10. Volumetric Rates 

An additional important factor for understanding the 
impact of potential contaminants on fast cadence, wide- 
held searches is their volumetric rates. We note that 


TABLE 2 

Expected Rates For Type Ia SNe 


Redshift 

(Mpc -3 yr -1 ) 

T^area 

(deg -2 yr -1 ) 

A 

0 - 0.25 

3 x 10 -5 

3.5 

21 

0.25 - 0.50 

4 x 10 -5 

25 

153 

0.5 - 0.75 

5 x 10 -5 

65 

394 


Note. — Expected rates for Type la SN e in redshift bins of 
A z = 0.25 out to Zmax = 0.75. See § |2.10| Graur et al. (2014), 
and references therein. 


the rates for these types of objects are generally poorly 
understood due to the lack of strong observational con¬ 
straints; the notable exception are Type Ia SNe. The 
volumetric rates of potential contaminants are compiled 
in Table [l] No values are provided in the literature for 
the ELDD volumetric rate so we assign a fiducial value 
of 1% of the Type Ia SN rate, consistent with the other 
rapid transients in this study. The volumetric rate for 
Type Ia SNe is known to be a function of redshift, and 
we use an approximate average volumetric rate in red¬ 
shift bins of Az ^ 0.25 out to z max ^ 0.75. (Li et al. 
2012, Perrett et al. 2013, Graur et al. 2014 and refer¬ 
ences therein). These rates are listed in Table [2j 
To estimate the expected number of events observed 
during a GW optical follow-up search we compute the 
number of sources per square degree per year. This is 
done by determining the maximum redshift to which the 
object can be obse rved to a lim iting magnitude of z ~ 23 
mag (see Sections 2.1 and |4.1| ). This rate is then multi¬ 
plied by the search area (A ~ 100 deg 2 ) and duration to 
get the number of expected sources for our observations. 
We compute the duration by assuming that a kilonova 
will be detected within approximately one day of the GW 
trigger. Therefore, we do not consider events that begin 
more than +2 d into the search. However, we do con¬ 
sider sources that begin up to twice their characteristic 
timescale (^ 5 d for rapid transients and ^ 20 d for Type 
Ia SNe) prior to the GW trigger, as they may remain vis¬ 
ible when the search begins. The total duration is then 
7 d for rapid transients and 22 d for Type Ia SNe. The 
number of expected sources is given in Tables [l] and [2] 
Type Ia SNe dominate the number of expected events 
with N ~ 568 when summing across all redshift bins. Of 
the population of fast transients, only type .Ia and the 
Pan-STARRS fast transients are expected to have N > 1. 


Despite the low expected occurrence rate for other fast 
transients, we still consider these sources in this analysis 
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Fig. 3. — Time-color (left) and time-magnitude (right) slices of the time-color-magnitude phase-space at a selection of redshifts. Arrow 
markers denote upper limits on the rise time. Note that in all panels the population of kilonovae are kept at the fiducial distance of ~ 200 
Mpc while redshift corrections are applied to the contaminant population. Kilonovae are well separated from the contaminant population 
in both color and timescales. The delineation in magnitude is less apparent. 


due to the large uncertainty in their rates. 

3. THE COLOR-MAGNITUDE-TIMESCALE PHASE-SPACE 

In this section we assess the possibility that any of the 
potential contaminants discussed in Section [2] could be 
confused for a kilonova during a GW follow-up observa¬ 
tion. We consider the maximum detection distance for 
a NS-NS binary merger to be 200 Mpc (z « 0.05), while 
for the contaminating sources we utilize larger distances 
as appropriate. We establish a three-dimensional phase- 
space based on the peak z-band brightness, the i — z 
color at z-band peak, and the rise time, defined as the 
time it takes the transient to rise one magnitude to peak 
in £-band. We compute the phase-space coordinates at 
200 Mpc for all of the transients discussed in Section [2] 
We then examine all of the potential contaminants at z 
~ 0.2, 0.4, and 0.6. In cases where time resolved spectra 
were available, the spectra were appropriately redshift- 
corrected in both time and wavelength. Sources mod¬ 
eled with a blackbody function had redshift corrections 


applied to the effective temperature. The rise time-color 
slice of this phase-space is shown in the left column of 
Figure [3] while the rise time-magnitude slice is shown in 
the right column of Figure [ 3 ] 

Investigating the first panel in the rise time-color space, 
we can see that the nine kilonovae models are well sep¬ 
arated from the population of contaminating sources by 
virtue of their redder color. We also note that their 
rise times are generally shorter, with only the AIC and 
1.0 + 0.05 M© type .la light curves exhibiting comparable 
timescales. The region occupied by kilonovae in the rise 
time-peak brightness phase-space is not as tightly con¬ 
strained. They are still well separated by timescale but 
several of the contaminating sources have comparable 
peak brightnesses at 200 Mpc. As we consider more dis¬ 
tant contaminant populations, we observe little to mod¬ 
erate reddening, with all of the contaminants remaining 
blue compared to the kilonovae population (i — z < 0 
mag to z ~ 0.6). The 0.6 + 0.3 M© and 1.2 + 0.05 M© 
Type .la SNe models have i — z ~ 0.2 mag, but their 
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longer timescales (t ~ week) keep them well separated 
from kilonovae. In addition, there is a systematic change 
in timescales, with all of the contaminants (except for 
the AIC events) shifting to a rise time of > 2 d by z ~ 
0.6 due to time dilation. 

The inherent low luminosity of most fast transients 
means that they quickly drop below the sensitivity limit 
of a reasonable search (e.g., m z « 23 AB mag). Only 
the Type la SNe, Pan-STARRS fast transients, and the 
most luminous type .la models are brighter than « 23 
AB mag by z ~ 0.2, and by z ~ 0.6 only the transient 
PSl-llqr and Type la SNe will still be detectable. 

Quantitatively, we can impose cuts of t r i se < 5 d, 
i — z > 0 mag, and m z < 23 mag to isolate kilonovae 
from contaminants. For sources at z < 0.2, this will 
isolate kilonovae from all contaminants except the Pan- 
STARRS fast transient PSl-lOah. These cuts remain 
effective for more distant sources as well. For example, 
at z ~ 0.6, AIC, the 0.6+0.3 M© Type .la SNe, and PSl- 
lOah all exhibit t r i se < 5 d and i — z > 0 mag. However, 
they no longer pass the brightness cut, with m z > 23 
mag. Therefore, we conclude that given well-sampled 
light curves, we can use the rise time-color-brightness 
phase-space to separate kilonovae from the contaminant 
population of other rapidly evolving sources. 

4. SIMULATED FOLLOW-UP OBSERVATIONS 

The analysis in the previous section was based on the 
idea that it is possible to obtain well-sampled light curves 
of kilonovae and other fast transients, but in reality the 
search cadence for a large area and the required depth 
will lead to at most a few light curve points. To ex¬ 
plore the impact of realistic observations we utilize Monte 
Carlo (MC) simulations which model a set of GW follow¬ 
up observations parameterized by cadence, limiting mag¬ 
nitude, and start time following a GW trigger. The 
goal of these simulations is two-fold. First, we study 
the feasibility of extracting useful light curve parame¬ 
ters for the kilonovae. This includes the actual detec¬ 
tion of a kilonova, as well as its detection on the rise, 
which is important for establishing the ejecta parame¬ 
ters. Second, we study the impact of the contaminant 
population on our ability to identify kilonovae by utiliz¬ 
ing our simulated observations to explore the rise time- 
color-magnitude phase-space established in Section [3j 

The detection of transient sources in wide fields ne¬ 
cessitates difference imaging. In the simplest scenario a 
pre-existing image of the field can be used as a template 
to eliminate the non-varying sources. In the context of 
real-time GW follow-up, however, such template images 
may not be available until after the source fades away. 
To address this difficulty, we consider three scenarios: 

1. The GW event is localized to a region with pre¬ 
existing template images. For example, the event 
occurs in a region previously observed to a similar 
or greater depth by the wide-held telescope being 
used for the follow-up. 

2. Template images are obtained after the source has 
faded. For kilonovae follow-up, this requires obser¬ 
vations at > 10 days after the GW trigger. 

3. The first observation taken as part of a GW follow¬ 
up campaign is used as the template image. 


From an analysis standpoint, the first two scenarios are 
identical as both involve a template image which contains 
no flux from the transient source. This allows robust 
subtractions to be performed, as well as precise mea¬ 
surements of brightness and color. However, the second 
scenario does not allow for the real-time identification 
of the transient, which is critical for obtaining follow-up 
spectroscopy and NIR imaging. These scenarios are the 
primary focus of this section. 

The third scenario, using the first observation as a tem¬ 
plate, is more challenging. The potential for the presence 
of flux from the transient in the template image means 
that we can only measure a relative difference in bright¬ 
ness. When this difference is small compared to the noise 
in the image it becomes difficult to measure key quanti¬ 
ties such as color or even claim that the source has been 
detected in difference imaging. We note that the second 
and third scenarios are not mutually exclusive. If the 
kilonova can not be identified in real-time under the third 
scenario, then a template image can still be obtained af¬ 
ter the source has faded away and scenario 2 becomes 
applicable. These issues are addressed in Section [5) 


4.1. Detectability of Kilonovae in Simulated 
Observations 


We first investigate the issue of kilonova detectability, 
as well as observation of the light curve during the rise 
to peak brightness. We explore four observing strategies 
utilizing ~1 week of observations in both i- and z-band 
with the following cadences: 

Strategy 1: One visit per night for the duration of the 
search (i.e., Nl, N2, N3, N4, N5, N6). 

Strategy 2: One visit per night on the first two nights, 
followed by a visit every other night (i.e., 
Nl, N2, N4, N6). 

Strategy 3: Two visits ~ 3 hr apart on the first night, 
followed by one visit per night for the du¬ 
ration of the search (i.e., Nl, Nl, N2, N3, 
N4, N5, N6). 

Strategy 4: Two visits ~ 3 hr apart on the first night, 
followed by a visit every other night (i.e., 
Nl, Nl, N3, N5) 


These strategies are chosen to probe a range of timescales 
while representing differing levels of time investment, and 
possible loss of coverage due to adverse weather during 
follow-up observations. 

For each observing strategy we establish a grid of start 
times following a GW trigger and a 5 a limiting mag¬ 
nitude. The range of grid values is £ sta rt = 3 — 24 hr 
in steps of 3 hr and mu m = 21 — 27 AB mag in steps 
of 1 AB mag. The range of start times represents the 
fact that a GW trigger may occur when the source is 
not readily visible in the night sky, while the range of 
limiting magnitudes covers the depths achievable by cur¬ 
rent and future wide-held instruments. At each point 
on this grid we conduct 50,000 simulated observations, 
with an observation being comprised of: (1) a kilonova 
with M e j = 10“ 2 M© and /3 e j = 0.2 placed at a distance 
of up to 200 Mpc using a uniform volume distribution; 
(2) the light curve of the kilonova is computed for the 
chos en d istance using the methodology outlined in Sec¬ 
tion 2. 1[ (3) the observation times are computed based 
on the chosen strategy; for example if the start time is 
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Fig. 4. — Depth versus start time grid in r (top), i (middle), and z (bottom) for the four observing strategies. The contours indicate the 
fraction of sources detected in at least three epochs (25% blue, 50% cyan, 75% yellow, 95% red). The dashed lines represent the approximate 
5 <t limiting magnitudes of several wide-held telescopes assuming the appropriate number of pointings necessary to cover a ~ 100 deg 2 error 
region with the required cadence. The choice of observing strategy does not drastically affect the fraction of sources detected at a given 
depth. We note that the early time drop in efficiency for z-band utilizing Strategy 2 is the result of the initial observations occurring before 
the light curve has brightened above the limiting magnitude. 


TABLE 3 

Required 5cr Limiting Magnitudes for Kilonova Detection 


Filter 

Strategy #1 

With Pre-Existing Template 
Strategy #2 Strategy #3 

Strategy #4 

Strategy #1 

Without Pre-Existing Template 
Strategy #2 Strategy #3 

Strategy #4 

r 24.9 (25.9) 

i 23.2 (23.9) 

z 22.2 (22.9) 

25.0 (25.9) 

23.3 (24.0) 

22.4 (23.2) 

24.5 (25.5) 
22.9 (23.9) 
22.1 (22.9) 

24.5 (25.5) 
22.9 (23.9) 
22.1 (23.0) 

25.0 (25.9) 
24.0 (24.9) 
23.1 (23.9) 

25.0 (25.9) 
24.3 (25.0) 
23.2 (24.0) 

25.0 (25.9) 
24.0 (24.9) 

23.1 (23.9) 

25.1 (25.9) 
24.3 (25.0) 
23.5 (24.2) 


Note. — Limiting magnitudes (5cr) required to detect 50% (95%) of kilonovae as a function of observing strategy and filter (Section |4.1| 
Figure pi. Values are averaged over all start times. We compute the required depths for follow-up observations with and without a 
pre-exisTing template (Section [5}. 


12 hours and strategy 1 is chosen, then observations will 
be conducted at 0.5 d, 1.5 d, 2.5 d, and so on; and (4) 
the source flux in rzz-bands DECam filters is measured 
at the selected times using the light curves computed in 
step (2). If the kilonova brightness exceeds the search 
depth in at least three epochs, the event is flagged as a 
detection. If there is an observed increase in brightness 
between two epochs, the event is also flagged as “rising.” 

In Figure [4] we plot contours of the fraction of events 
detected for each of the four observing strategies in the 
riz- bands, as a function of survey depth and start time. 
The dashed lines represent the approximate 5 a limiting 
magnitudes of several wide-held telescopes assuming the 
appropriate number of pointings necessary to cover a ^ 
100 deg 2 error region with the required cadence. 

We find that the choice of observing strategy does 
not strongly impact the fraction of sources detected at 
a given depth. In particular, in r-band (Figure [4j top 
row) we find that for both Strategies 1 and 2, the typical 


depth required to detect 50% of events is ~ 25 mag, while 
a depth of ~ 26 mag is required to detect 95% of events. 
This challenge is slightly alleviated by Strategies 3 and 
4, which involve two observations on the first night. In 
this case, an early start time (< 6 hr) allows 50% (95%) 
of events to be detected at a depth of ~ 24.3 (25) mag. 
Consequently only Subaru HSC and LSST will be able 
to achieve the required depth for routine detections in 
r-band. The redder bands offer more favorable condi¬ 
tions, with the z-band 50% (95%) detection rate achieved 
at a depth of « 23(24) mag across all four strategies. 
Lastly, we find that the z-band contours indicate that a 
depth of ~ 22.2 (23) mag is required to achieve a 50% 
(95%) detection rate across all four strategies. In these 
cases, DECam, Pan-STARRS2, and Subaru HSC can all 
achieve the required depths for routine detections in i- 
and z-bands. These values are summarized in Table 02 
We also find that for z-band the observing strategies 
utilizing longer cadences (i.e., Strategies 2 and 4) produce 
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Fig. 5. — Same as Figure^] except that the contours show the fraction of events for which a brightening is detected. The choice of 
observing strategy, start time, and filter have a strong impact on the detection of brightening. 


a noticeable drop in efficiency if the first observation oc¬ 
curs too early (< 6 hr). In this case, the light curve 
has not yet brightened above the survey depth, causing 
the initial observation to be a non-detection. The slow 
cadence after the initial night of observations (at > 2 
days) is then insufficient to capture the requisite three 
detections before the light curve fades below the detec¬ 
tion limit. Quantitatively, we find a loss of ^ 1 mag for 
Strategies 2 and 4 for starting times < 6 hr. This indi¬ 
cates that the strategies involving slower cadence after 
the initial night of observations are not well matched to 
kilonova timescales. 

While the particular choice of search strategy does not 
strongly impact the detection fraction at a given survey 
depth, the choice of strategy is important when trying 
to detect the brightening of a kilonova. We investigate 
this point in Figure [5j which shows contours for the frac¬ 
tion of events that are seen to brighten between any two 
epochs. These events must also meet the detection crite¬ 
rion utilized in Figure [4] As in Figure [4j the contours are 
plotted as a function of start time and search depth for 
each of the observing strategies and filters. In r-band, 
utilizing Strategies 1 and 2 , a brightening phase is not 
detected even to a limiting magnitude of 27 mag. This 
is because the r-band light curve peaks in < 1 day (Fig¬ 
ure [l]) necessitating a rapid early time cadence to capture 
the event as it brightens. For Strategies 3 and 4, with 
two observations on the first night, brightening can be 
observed for 50% (95%) of events provided the observa¬ 
tions reach a limiting magnitude of ~ 24.3 (25) mag and 
begin within 8 ( 6 ) hr of the GW trigger. 

The detection of the light curve rise is more favorable 
in i- and £-band. In general, the contours shown in Fig¬ 
ure [5] are identical to those in Figure [4j but truncated as 
the start time becomes comparable to the peak time of 


the light curve. Therefore, brightening can be detected 
provided observations begin early enough and assuming 
the required depth is achieved. For z-band, with Strate¬ 
gies 1 and 2 , we find that observations starting within 6 
( 8 ) hr of the GW trigger can detect brightening in 50% 
(95%) of events. With Strategies 3 and 4 the required 
start time for a detection rate of 50% (95%) extends out 
to 12 (14) hr. In z-band the observations can achieve 
detection rates of 50% (95%) provided they begin within 
12 (15) hr of the GW trigger. The starting time can be 
extended to 18 ( 21 ) hr when Strategies 3 and 4 are uti¬ 
lized. We also note that Strategies 2 and 4 show the same 
early time (< 6 hr) loss of ^ 1 mag as seen in Figure [4j 
The broader range of starting time allowed by strategies 
3 and 4 highlight the advantage of rapid follow-up of GW 
triggers. However, these initial simulations consider any 
observed brightening as the detection of a rise. It may 
be the case that the actual change in brightness is not 
statistically significant, making a robust detection of a 
rise impossible. 

We investigate the actual change in brightness rela¬ 
tive to the instrument sensitivity using a modified ver¬ 
sion of our initial simulations. We no longer operate on 
the depth-start time grid, but instead consider a fidu¬ 
cial fixed depth of 24 AB mag in r- and z-band and 
23 AB mag in z-band (as appropriate for achieving an 
acceptable detection fraction). The start time is ran¬ 
domly selected using a uniform distribution ranging from 
3 to 24 hr. The same four observing strategies are em¬ 
ployed and we run 50,000 iterations per strategy, uti¬ 
lizing the same procedure as before. We then compute 
the change in magnitude between the first two epochs 
(A? 7 Z 2 i = m 2 — mi), and between the first and third 
epochs (Am 3 i = m 3 — mi). We employ the same detec¬ 
tion criteria as in the previous simulations, requiring the 
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Fig. 6. — Probability density of kilonovae as a function of the change in brightness between the first and second observation (A 7712 I) versus 
the brightness in the second observation (m 2 ) for each of the four observing strategies. The hatched region indicates the 2cf photometric 
error for DECam for an exposure time necessary to achieve the required depth. Points within this region will not exhibit a statistically 
significant change in magnitude for detection. The fraction of sources brighter than the detection limit in all three observations (DF) and 
the fraction of sources within the hatched region (missed) are indicated at the bottom of each panel. We find that the seemingly beneficial 
initial rapid cadence of Strategies 3 and 4 on a ~3 hr baseline actually result in a large fraction of brightenings being undetectable. 


source to be detected in at least three epochs. We con¬ 
sider sources that are observed to either rise or decline. 

In Figure [ 6 ] we plot two-dimensional histograms of 
A 777/21 versus m 2 for each of our four observing strategies 
and each of the three filters. The resulting histogram 
is color-coded by the probability of a kilonova falling in 
that bin, with darker colors indicating a higher proba¬ 
bility. We also compute the 2 cr photometric error for 
an exposure time necessary to reach our chosen target 
depth in each filter using the DECam exposure calcu- 
latoi 0 (hatched region). Sources that exhibit a change 
in magnitude that is smaller than this photometric error 
will not produce a statistically significant detection of a 
change in brightness. Figure [7] follows a similar construc¬ 
tion, showing the two-dimensional histograms for A 777/31 
versus m 3 . 

Using these results we determine the fraction of sources 
that do not exhibit a sufficient change in brightness to 
be identified as a transient after the second or third ob¬ 
servation. We begin by considering the first and second 
columns of Figures [ 6 ] and [7] (Strategies 1 and 2 ). Recall 
that for these strategies the time between the first and 
second epoch is ~1 day. In both of these cases, the rapid 
evolution of the r-band light curve means that there is 
always a detectable change in brightness, but the source 
will always be seen in decline. Furthermore, it is still 
the case that given the overall faintness of kilonovae in 
r-band, only ^13% of sources will be detected at r ~ 24 
mag. 

In i- and z -band we find that (^ 20 — 40%) of sources 
are missed due to an insufficient change in brightness. 
For Strategies 3 and 4 where the time between the first 

;) http://www.ctio.noao.edu/noao/content/ 

Exposure-Time-Calculator-ETC-0 


two observations is only ^3 hr over half of the sources 
(60 — 70%) fail to exhibit a statistically detectable change 
in brightness. This indicates that the intra-night obser¬ 
vations have a time baseline that is too short for de¬ 
tecting a change in brightness even for the fast kilonova 
timescale. 

Figure [7] shows observations that are separated by a 
longer timescale (m 3 — mi): For Strategies 1 and 4, the 
first and third epoch are separated by 2 days, while the 
same visits are separated by 4 days for Strategy 2 , and 
only 1 day for Strategy 3. For all four strategies we 
find that r-band exhibits a sufficient change in bright¬ 
ness, with Am 3 i well outside of the hatched region, but 
these sources are always observed to decline. As before, 
only ^ 10 — 20 % of the sources are actually detectable at 
r ~ 24 mag. In both i- and £-band we find that it is pos¬ 
sible to detect a change of brightness in a much higher 
fraction of sources compared to Am 2 i- Strategies 1 and 
4 both show that ^ 15 — 20% of sources will not show 
a sufficient change in brightness. The highest fraction of 
missed sources is ^ 25 — 40% for Strategy 3. This is ex¬ 
pected given the short separation between epochs. Strat¬ 
egy 2 , which has the longest separation between epochs, 
shows only ^ 10—15% of sources remaining undetectable. 
Furthermore, while sources in r-band are only observed 
to decline, in i- and 2 -band the rise can also be mea¬ 
sured without requiring a ^3 hr baseline. This further 
highlights the importance of having observations that are 
sufficiently separated in time to probe the kilonova light 
curve across a wide range of its evolution. 

The rapid decline of ~ 0.6 — 1.3 mag in r-band that is 
evident in Figures [ 6 ] and [7] hints at a possible detection 
strategy. Namely, the GW error region can be observed 
with a cadence of d, with rapidly declining sources 
flagged. However, a depth of ~26 mag is required for a 

































































Contamination in GW Optical Follow-up 


11 


Strategy #1 


-0.4 

- 0.8 


Strategy #2 



Strategy #3 


Strategy #4 





'■r - 


r 

PF: 19 %, Missed; 0% 


PF: 13 %, Missed; 0 % 


0.8 
0.6 
0.4 
0.2 
0.0 
- 0.2 
-0.4 

-0.6 - DF: 100 %, Missed: 14 %j [, DF: 98 %, Missed: 9 °/<s\ ^ DF: 100 %, Missed: 24 %\ ^ DF: 100 %, Missed: 14 °/<r 

0.8 


r .. 



r 




I 


24 22 20 18 16 14 24 22 20 18 16 14 24 22 20 18 16 14 24 22 20 18 16 14 

m 3 [Mag] m 3 [Mag] m 3 [Mag] m 3 [Mag] 


> 

!75 

c 

CD 

o 

>> 

fO 

_Q 

O 

CL 


Fig. 7.— Same as Figure [ 6 ] except we plot the change in brightness between the first and third epoch (Am 3 i) versus the brightness in 
the third epoch (m 3 ). We find that the increased time between observations allows a much larger fraction of sources to exhibit a detectable 
change in brightness. 


high detection fraction even if we require only 2 detec¬ 
tions (to separate kilonovae from stellar flares), so cover¬ 
ing the large GW error region can only be achieved with 
LSST. 

These detectability simulations indicate that using 
intra-night observations (Strategies 3 and 4) to detect 
the rise of the light curve is not an efficient use of time 
for follow-up observations since kilonovae generally do 
not exhibit a statistically significant change in bright¬ 
ness on this timescale. As shown in Figure [4j the choice 
of cadence does not drastically impact the likelihood of 
detection. Therefore, the most efficient solution, from a 
detectability standpoint, is Strategy 1. It provides a suf¬ 
ficient cadence to effectively detect kilonovae while max¬ 
imizing the amount of area that can be covered in a sin¬ 
gle night. Observations should be carried out in i- and 
z-band as the depths required to detect 95% of sources 
(~ 24,23 AB mag, respectively), are obtainable by ex¬ 
isting and future facilities. 

4.2. Simulating the Effects of Contamination 

We now use our simulated observations to explore the 
contaminant population in the framework of the rise 
time-color-magnitude phase-space outlined in Section [3j 
We accomplish this by repeating the simulations for each 
of the contaminants considered in Section [21 We then uti¬ 
lize the simulated observations to compute analogues to 
the rise time-color-magnitude phase-space for both kilo¬ 
novae and the contaminant population. In this manner 
we can produce a version of the phase-space that takes 
into account a distribution of contaminants with realis¬ 
tically sampled light curves given the choice of cadence, 
depth, and start time appropriate to wide-held follow-up 
observations of a GW trigger. 

We begin by employing the same four observing strate¬ 
gies as in the previous simulations. We again consider 
observations at a fixed depth of ~ 24 AB mag in z-band 


and « 23 AB mag in z-band, with the observations be¬ 
ginning 3 to 24 hr after the GW trigger. The simulations 
include 50,000 iterations, each consisting of: (1) A con¬ 
taminant placed at a random, volume-weighted redshift. 
We consider contaminants at much larger distances than 
200 Mpc to account for the fact that many contaminants 
are more luminous than kilonovae; the choice of z max is 
indicated in Figure [8j (2) The contaminant light curve 
is randomly placed out of phase with respect to the start 
time of observations. We consider contaminants that 
occur approximately twice their characteristic timescale 
prior (5-20 days) and up to 2 days after the start of ob¬ 
servations. The latter value is motivated by the results 
of the previous section in which we showed that kilono¬ 
vae should be detectable within 1-2 days, and therefore, 
sources that first appear more than 2 days after the start 
of the observations can be rejected. Once the phase off¬ 
set is determined the observation times are computed 
and the source brightness in i- and z-band is computed 
at each epoch. (3) The observed peak of the z-band light 
curve is measured. (4) The i - z color is computed at the 
epoch where z-band attains its maximum value. (5) The 
rise time is defined as the time between the z-band peak 
and the first observation in which the source is detected. 
In cases where the source is only observed in decline, then 
the decline time is defined as the time between the first 
observation and the last observation in which the source 
is detected. 

In Figure [5] we plot the rise time-color slice of the three- 
dimensional phase-space for each observing strategy. In 
addition to the nine kilonovae models, we consider four 
groups of sources with each given an equal contribution. 
The first group of contaminants is Type la SNe (Figure[8j 
first row). The second group is comprised of the four type 
.la SNe models (Figure|8j second row). We form the third 
group from all of the low luminosity, fast transients (i.e., 
AIC, WD-NS/BH mergers, and ELDD; Figure [8| third 
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Fig. 8. — Simulated rise time-color phase-space for each of our four observing strategies. In each panel, we consider the full grid of 
kilonova models to a distance of 200 Mpc (orange). Each row highlights a different subset of the contaminant population over the comoving 
volume to z max (blue). The offset bins for Strategies 3 and 4 are due to the availability of intra-night observations on the first night. We 
find that the kilonovae are much redder than any potential contaminant, but the separation is less well-defined when compared to the light 
curves used in Section [3] 


row). Lastly, the fourth group of contaminants contains 
the three Pan-STARRS fast transients (Figure [8j fourth 
row). 

To analyze the results, we first bin sources by their rise 
times, which fall in discrete values due to the cadence of 
the observations. For each rise time bin, we then bin the 
source populations by i — z color in equally spaced bins 
from -1.5 to +1.5 mag. For each bin, we compute the 
fraction of detected sources in that bin relative to the 
total number of detected sources. 

As we showed in Section [3j kilonovae stand out from 
the contaminant population primarily on the basis of 
their redder colors, even when luminous contaminants 
are considered to higher redshifts (e.g., Type la SNe, 
Pan-STARRS fast transients). For all choices of observ¬ 
ing strategy there is essentially no overlap in i — z color 
between the kilonova population and any of the contami¬ 
nant populations. We note that across all four observing 
strategies ^12% of Type la SNe, ^ 7% of type .la, and 
~ 8% of Pan-STARRS fast transients exhibit i — z > 0, 
while none have i — z > 0.4 mag. However, all kilonovae 
exhibit i — z > 0.4 mag. These numbers are relative to 
the total number of events observed during the rise and 
do not include sources that are observed only in decline 


(see below). Consequently, i — z color remains a strong 
discriminant for separating kilonovae from the contami¬ 
nant population. 

So far we have neglected the potential effect of extinc¬ 
tion in reddening the i — z colors of the contaminant 
sources. If the extinction is sufficiently large an intrinsi¬ 
cally blue contaminant may be reddened into the regime 
of kilonova i — z colors. However, this extinction will 
also cause the source to appear dimmer and may there¬ 
fore push it below the detection limit. To quantify the ef¬ 
fects of extinction we use the Milky Way extinction curve 
with Ry =3.1. Figure [8] indicates that a reddening of 
E(i — z) ~ 1 mag is required to shift the contaminant 
population into the region of phase-space occupied by 
kilonovae. This will result in an extinction of Ay ^6.3 
mag, Ai « 3.9 mag and A z « 2.9 mag. However, in our 
simulation nearly all (> 97%) of the contaminants which 
pass the detection cuts have m z > 20 mag meaning that 
an extinction of 2.9 mag will dim these sources below our 
z-band detection limit of ~ 23 mag. Thus, reddening due 
to dust extinction does not seem to be a concern. 

We also consider the possibility of separating kilonovae 
by their rise time. In Figure [3] the entire set of kilono¬ 
vae models showed generally shorter timescales relative 
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Fig. 9. — Same as Figure[8] but showing the simulated rise time-A m z slice of the phase-space. In this slice, there is much more confusion 
between the kilonovae and contaminant populations. With the exception of Type la SNe, there is significant overlap between kilonovae 
and the contaminant populations. This is the result of considering rapidly evolving sources for our contaminants, and it highlights the 
importance of color information in the selection of kilonovae candidates. 


to the contaminant population with rise times of < 4 
days. However, when the contaminants are placed out 
of phase with the GW trigger, and the observing ca¬ 
dence is taken into account, we find contaminants that 
present an apparent shorter rise time. This effect is the 
most prominent in the Pan-STARRS fast transients and 
miscellaneous fast transients for which ~ 55 — 60% and 
^ 65 — 75% of detected sources show a rise time of < 4 
days, respectively. We also find that ^ 35 — 45% of de¬ 
tected type .la show a rise time that is < 4 days. Lastly, 
for Type la SNe, ^ 20 — 25% of detected sources will 
present a rise time of < 4 days. Thus, rise time alone is 
not a strong discriminant for separating kilonovae from 
the contaminant population. 

In Section [3J we considered the peak z -band magni¬ 
tude as the third data point in our phase-space. How¬ 
ever, in a magnitude limited search, the population of 
contaminants will be dominated by events at the fur¬ 
thest distance to which they can be detected. Here we 
instead track the evolution of the source brightness by 
computing the difference in z -band magnitude between 
the source at the observed peak (m^ jPea k) and the first 
observation in which it is detected (m Zj o), denoted as 
A m z = ra^peak — m z ^. We show the rise time-A m z slice 
of the three-dimensional phase-space in Figure [9j The 


basic layout of the figure is identical to that of Figure [8j 
but sources are offset by ±5 hours to improve visibility. 

We find that kilonovae exhibit A m z « —3.5 to 0 mag, 
where larger changes in brightness are associated with 
longer rise times. Type la SNe only exhibit values of 
A m z comparable to those of kilonova when the observed 
rise time is ^ 6 d (the duration of our search), well 
beyond the range of observed rise times for kilonovae. 
When the rise time is < 6 d, the average change in bright¬ 
ness is negligible (A m z « 0.05 mag). This is because 
shorter rise times only occur near peak when the light 
curve evolution is slow. This slow evolution will make 
Type la SNe easy to reject relative to the much more 
rapidly evolving kilonovae. We note a small fraction of 
events with A m z « — 2 mag and rise times of 4 — 5 days. 
These account for < 0.2% of detected Type la SNe and 
are the result of an event occurring in a range of redshifts 
where the set of spectral features that produce the sec¬ 
ond peak of Type la SNe appear in z-band. Thus, Type 
la SNe can be rejected based on A m z , rise time, and i — z 
color. For the other contaminant populations the sepa¬ 
ration is not as clean. This is not unexpected since rapid 
evolution was a strong consideration when constructing 
our list of potential contaminants. This highlights the 
critical importance of selecting kilonovae on the basis of 
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Fig. 10. — The decline time-color slice of the decline time-color-magnitude phase-space. The plot is constructed using the same method¬ 
ology as in Figures [8] and [9] We note that, as before, kilonovae show redder i — z colors than the contaminant population. However, there 
is no strong separation in aecline time. 


their i — z color along with their rapid rise times. 

We demonstrated in Figure [5] that it is challenging to 
detect a kilonova during the rise to peak. Therefore, in 
our simulations we also consider the behavior of sources 
that are observed in decline. We construct an identical 
phase-space by replacing the rise time with the decline 
time, as defined above. Figure [10] shows the decline time- 
color slice of this phase-space. The layout is identical to 
Figure [ 8 | As before, the kilonova population presents 
i — z colors that are consistently redder than those of 
the contaminant population with all detected events ex¬ 
hibiting i — z > 0.3 mag. Furthermore, the fraction of 
contaminants seen in decline that exhibit i — z > 0 is es¬ 
sentially unchanged compared to those observed on the 
rise. Additionally, all of the contaminants considered still 
exhibit i — z < 0.3 mag, so therefore color can still be 
utilized as a strong discriminant. 

The efficacy of a cut on decline time is diminished com¬ 
pared to a cut on rise time. Both kilonovae and the 
contaminant population cover the entire distribution of 
possible decline times from ^2 — 6 days. A notable ex¬ 
ception is Type la SNe, which all exhibit a decline time 
of ^ 6 days (the duration of the search), because if the 
source is already in decline at the start of the observa¬ 
tions, it will remain detectable for the entire duration of 
the search. Only ^ 20—25% of kilonovae exhibit a decline 


time of < 6 days, making this an ineffectual cut. On the 
other hand, the slow evolution is beneficial for separating 
the Type la SNe population from a kilonova if template 
images at > 10 days are utilized. In this scenario, the 
kilonova will have faded away while Type la SNe will 
remain visible allowing them to be easily rejected. 

Figure [Tl] shows the decline time- A m z slice of the 
phase-space. In this context we define A m z = m Zj i ast — 
m Zj o, where m Z: i ast is the £-band magnitude from the last 
observation in which the source is detected. Kilonovae 
exhibit a change in brightness of A m z ^0 — 4 mag, with 
a weak trend for larger changes in brightness being as¬ 
sociated with longer observed decline times. The Type 
la SNe and Pan-STARRS fast transients exhibit an av¬ 
erage change in brightness of only A m z « 0.1 mag and 
A m z « 0.2 mag, respectively, making the potential for 
contamination low. However, as with the rise time-A m z 
space, the type .la and miscellaneous fast transients ex¬ 
hibit changes in brightness that are comparable to those 
of kilonovae. This makes cuts on A m z less practical for 
these sources. 

The key point emerging from our simulations is that 
i — z color is a strong discriminant for selecting kilonovae, 
regardless of observing strategy, and whether a kilonova 
is observed on the rise or decline. In fact, if we consider 
sources observed during the rise to peak, and we require 
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Fig. 11. — The decline time-A m z slice of the decline time-color-magnitude phase-space. The slow decline of Type la SNe and the Pan- 
STARRS fast transients separate them from the rapidly evolving kilonovae. However, the type .la and the miscellaneous fast transients 
show changes in brightness comparable to those of kilonovae. 


that i — z > 0 mag we may eliminate ^ 84 — 87% of 
detected Type la SNe and Pan-STARRS fast transients, 
^ 90 — 94% of detected type .la, and ^ 100% of detected 
miscellaneous fast transients without removing any kilo¬ 
novae from the selection. If we additionally require that 
the rise time is < 4 days, then we eliminate ^ 94% of 
detected Type la SNe and ^ 98% of detected type .la. 
If we also require that lAraJ > 0.1 mag, we eliminate 
- 99% of Type la SNe. 

Considering sources observed in decline the cuts on i—z 
color and A m z remain effective but a cut on decline time 
does not provide any additional information. Quantita¬ 
tively, if we apply cuts of i — z > 0 mag and | A m z | > 0.1 
mag we eliminate ^ 80% of Type la SNe, ^ 60% of type 
.la, ~ 81% of Pan-STARRS fast transients, and ~ 100 % 
of the miscellaneous fast transients. If we also require 
that the decline time is < 5 days we eliminate ^ 100% 
of the Type la SNe and ~ 94% of the Pan-STARRS fast 
transients, but lose ^ 46% of the kilonova population. 
However, if additional observations are taken at > 10 
days then these kilonovae can be recovered. But, given 
that cuts made on sources observed to rise are more effec¬ 
tive, then in order to maximize the chances of a kilonova 
detection and identification during a GW follow-up cam¬ 
paign it is essential to obtain both i- and z -band data as 
close to the GW trigger as possible. 


5. COUNTERPART IDENTIFICATION IN THE ABSENCE 
OF A TEMPLATE 

We now investigate the situation where there are no 
pre-existing template images of the GW error region and 
we would like to rapidly identify a kilonovae without 
having to wait for a late-time template. This scenario 
introduces several challenges. Foremost is the presence 
of source flux in the image that will be used as a tem¬ 
plate, which will have an effect on the observables (e.g., 
brightness, color). Specifically, we can no longer mea¬ 
sure the true flux of the source but rather only the dif¬ 
ference in flux between the science and template images. 
In this case, a measurement of the true magnitude of the 
source is impossible. Without the accurate measurement 
of source brightness it is also difficult to measure the true 
color and we therefore lose one of the primary methods 
of distinguishing kilonovae from the contaminant pop¬ 
ulation. As discussed in detail below, these difficulties 
can be alleviated by making simple assumptions about 
the source behavior. 

Another obvious challenge is a delay in the timescale 
of the first science image since the first observation must 
be used as the template image. This is leads to problems 
in characterizing the early behavior of the light curve, 
as well as in utilizing the rise time as a method of dis¬ 
tinguishing kilonvoae from the contaminant population. 
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Fig. 12 .— Contours showing the fraction of kilonovae detected as a function of start time and depth for our four observing strategies and 
three filters when a pre-existing template is not available. We first note that i- and z-band observations in this scenario must be ~ 1 mag 
deeper to achieve detection rates comparable to when a pre-existing template is available. We also note an effect in z-band where the 95% 
contour drops by ~ 1 mag for observations starting between 6 and 12 hours after trigger using Strategy 2 and for observations starting 
less than 12 hours after the trigger using Strategy 4. This is the result of the observations straddling the peak of the light curve where the 
difference in flux between epochs is too small to detect. 


Using the first observation as a template also reduces 
the number of available epochs making it more difficult 
to meet the detection criteria outlined in the previous 
section. However, these challenges are less severe as they 
can resolved using a late-time template once potential 
counterparts are identified. 

To address the scenario of using the first image as a 
template, and to gauge the associated complications, we 
simulate obser vatio ns using the same methodology out¬ 
lined in Section |4T| but taking into account the fact that 
the first observation is the template. The general pro¬ 
cedure is identical to the previous simulations with the 
following modifications: We now simulate the effects of 
difference imaging by subtracting the flux in the first ob¬ 
servation (the template) from the later observations (the 
science images) to produce a set of difference fluxes. We 
account for the noise appropriately by using our search 
depth as a 5cr flux limit. We assume that this noise is the 
dominant source of error and that it is the same for all 
images. The noise in the difference image is then larger 
by a factor of y/2 than that of the individual images. 
Any source that shows a change in flux with a signal-to- 
noise ratio (SNR) greater than 5 in the difference image 
in at least two subtractions is flagged as a detection. If a 
source exhibits at least one epoch with a positive change 
in flux, the event is flagged as “rising.” 

In Figure [l2j we plot contours of the fraction of sources 
that pass our detection criterion as a function of start 
time and search depth. The layout of this plot is iden¬ 
tical to that of Figure [4] We first note that the r-band 
contours for Strategies 1 and 2 are essentially identical 
to those in Figure [4j This is because the r-band light 


curve evolves fast enough between observations that the 
effects of using the first observation as a template does 
not hamper detections. However, for Strategies 3 and 4, 
we no longer find the early-time (< 6 hours) boost in de¬ 
tection limit. This is due to the fact that the light curve 
does not change sufficiently during the three hours, and 
the source is therefore not detected in the first difference 
image. There is therefore no benefit from the additional 
observation on the first night. On the other hand, in 
z-band we find that on average the depth required to 
reach a certain detection rate is ~ 1 mag deeper than 
in the case of a pre-existing or late-time template (Fig¬ 
ure^. Specifically, for Strategy 1 the 50% contour starts 
at « 24.4 mag for £ start ^ 3 hours and rises to « 23.7 
mag by £ star t ^ 24 hours. In Strategy 2, the 50% contour 
is flat at « 24.3 mag. The 95% contour is flat at « 25 
mag in both cases. Strategy 3 is identical to Strategy 1, 
again highlighting the fact that the rapid cadence on the 
first night does not boost the detection rate. Lastly, we 
find that the rapid early-time cadence combined with the 
slow late-time cadence of Strategy 4 actually produces a 
drop in efficiency for observations beginning at ^ 3 hours 
where the 95% contour dips to ~ 25.7 mag. 

The most pronounced change relative to the case of 
a pre-existing template (Figure HI) can be seen in the 
band contours. In Strategy 1, the 50% (95%) contour 
is approximately « 23 (23.9) mag, which is again ~ 1 
mag deeper than for the case of a pre-existing template. 
Strategy 2 shows a change of behavior in the 95% contour 
between 6 and 12 hours, where the required depth dips 
to ~ 24.8 mag. This is because observations starting 
during this time frame will be roughly symmetrical about 
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Fig. 13.— Same as Figur e 1 12| but showing the fraction of sources for which brightening is detected. We note that when compared to the 
contours shown in Figure [5] for the case of a pre-existing template, observations must go deeper in this scenario and the detection efficiency 
quickly tapers off as a function of starting time. This highlights the importance of commencing follow-up observations as soon as possible 
when a pre-existing template is not available. 


the peak of the light curve. Therefore, the change in 
flux is small and the source will not be detected in the 
difference image. The slow cadence of Strategy 2 then 
makes it difficult to obtain enough observations to reach 
our threshold for a detection. This effect is also seen in 
Strategy 4 where the 95% contour dips to « 24.8 mag 
for any starting time less than 12 hours. These contours 
are summarized in Table 03 

In Figure [l3j we show contours for the fraction of 
sources with a detected brightening between the template 
and at least one of the science images. These sources 
must also pass the detection criterion described above. 
The layout is identical to that of Figure [5] The behav¬ 
ior of the contours is different compared to the case of a 
pre-existing template for all bandpasses and strategies. 
In Figure [5j the contours were flat and then exhibited 
a sharp drop off when the first observation began post 
peak. Here, the contours show a steady decline as a func¬ 
tion of start time. This is because as the start time of 
the first observation (i.e., the template) approaches the 
peak, the observed change in flux will decrease. There¬ 
fore, it becomes progressively more difficult to detect a 
brightening and the required depth increases. 

Specifically, for r-band, no brightening is detected 
when utilizing Strategies 1 and 2. The same is effec¬ 
tively true for Strategies 3 and 4, unless the observa¬ 
tions reach a depth of ~ 26.2 mag and begin < 3 hr 
after the GW trigger. Thus, even with deep Subaru 
HSC or LSST observations, brightening could only be de¬ 
tected in ~ 25 — 50% of sources for observations starting 
within 3 hours of the GW trigger. In z-band, for Strate¬ 
gies 1 and 2, the 50% (95%) contour begins at ~ 24.4 
(25.3) mag and no rise is detected if 4tart > 7 (6) hours. 
For Strategies 3 and 4, the 50% (95%) contours start at 


« 24.5 (25.3) mag and fall below the LSST limiting mag¬ 
nitude by ^ 9 (6) hours. Thus, the detection of a rise in 
z-band is not feasible with existing telescopes without a 
pre-existing template. Lastly, for z-band, the 50% (95%) 
contours in Strategies 1 and 2 begin at « 22.7 (23.9) mag 
and fall below the LSST limiting magnitude by ~ 15 (12) 
hours. For Strategy 3, the required depths are unchanged 
compared to Strategies 1 and 2, but the 50%(95%) con¬ 
tours extend out to ^ 17 (15) hours. For Strategy 4, the 
required depth to achieve a 50% (95%) detection rate 
is ~ 1 magnitude deeper compared to Strategy 3, but 
the required start times are identical. As above, this 
is possible with Subaru HSC and LSST, but other in¬ 
struments (e.g., DECam, Pan-STARRS2) will only be 
able to achieve a 95% detection rate for z-band observa¬ 
tions starting within ~ 6 hours of the GW trigger. This 
highlights the increased difficulty in characterizing the 
early-time behavior of the light curve because we lose 
information from the first epoch when it is used as the 
template image. We note that this information can be 
recovered after the fact using a late-time template but 
here we are concerned with rapid identification. 

5.1. Contaminant Rejection In The Absence of a 
Template 

We now investigate how the absence of a pre-existing 
template affects our ability to reject contaminants. We 
carry out our simulations in the following way. The 
flux in the first observation (the template) is subtracted 
off from the later observations. Noise calculations are 
handled identically to those outlined above for the de¬ 
tectability simulation, and we impose the same detection 
criterion for events. The rise/decline time phase-space 
coordinates are computed using the same definition as 
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Fig. 14. — The rise time-color slice of our phase-space when a pre-existing template is not available. We note that if the rise time is long 
(> 2 days), then there is sufficient evolution in the light curve that colors can be measured with reasonable accuracy and kilonovae can be 
separated from the contaminants as before. However, when the rise time is short (~ 1 day), then the change in flux is small and the colors 
are not well constrained. 


4.2 The challenge then is computing an ana- 


in Section 

logue for the values of the i — z color and A m z when the 
source flux is not known absolutely (due to the presence 
of source flux in the template image). 

When there is source flux in the template image we 
are only able to measure a difference in flux between 
the template and science image. We denote this as 
A /0 = / s ,0 - f t ,cf> where f s ^ is the flux observed in 
the science image and f t ,0 is the flux observed in the 
template image, both in some bandpass 0. We can 
then define an AB magnitude using the standard formula 
rat = — 2.5 log | A/ 0 I — 48.6, where the star denotes mag¬ 


nitudes computed from a difference in flux. The 0 1 — 02 
color for bandpasses 0 1 and 02 is then ra^ — ra^ 2 . In the 
case that | A/^ | is large, and consequently the detection 
is of high statistical significance, we approach the limit 
where the source flux in the template is negligible and 
a reasonable approximation for the true magnitude and 
color of the source can be obtained. However, for less 
significant detections, the measured color will not be an 
accurate measurement of its true value. 

Lastly, we construct a straight-forward analogue to the 
rise/decline time-Ara slice of our phase-space. In the 
case of a rising source, we compute raj using the largest 
positive difference in flux between the template image 


and the science images. If the source is observed only in 
decline then the largest negative difference in flux is used. 
This will allow us to understand how rapidly a source 
evolves relative to others in this study (e.g., kilonovae 
versus contaminants). 

In Figure [l4| we plot the rise time-color slice using the 
same layout as in Figure [8j Inspecting the kilonovae 
across all four observing strategies we note that when 
the observed rise time is small (< 1 day, i.e. when the 
difference in time between the template and first science 
image is small) then the color is not well constrained 
due to | A/1 being too small to obtain an accurate mea¬ 
surement. However, the fraction of kilonovae with a rise 
time of < 1 day is only ^ 5 — 8 %. We also note that 
the fraction of detected kilonovae with t r i se < 1 day that 
exhibit i — z ~ — 2 to +0.5 mag (i.e., in the range that 
overlaps with contaminants) is at most ^ 1 % for Strate¬ 
gies 3 and 4, and < 1% for Strategies 1 and 2 . More 
importantly, when the observed rise time is longer (> 2 
days), then |A/| is large enough that the i — z color can 
be well constrained and it is generally i — z > 1 mag. 
This is consistent with the color measured when a pre¬ 
existing template is available (Figure [ 8 ]). The fraction of 
detected kilonovae that present a rise time of > 2 days 
is ^ 45 — 55%. This fraction is dominated by kilonovae 
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Fig. 15. — The rise time-mj slice of our phase-space. The slower evolution of Type la SNe is apparent. For the remaining contaminants 
the observed change in flux is comparable to that observed for kilonovae. The rise times exhibited by both the kilonovae and contaminant 
populations are also comparable, however only contaminants show a rise time > 5 days. 


models with an ejecta mass of M e ,j = 10 _1 M 0 and /3 e j = 
0.1 or 0.2. 

Turning to the contaminant population, we find that 
several of the source classes have poorly constrained col¬ 
ors with some events exhibiting i — z > 0 mag. We can 
assess the impact of these contaminants by utilizing the 
rate data in Tables [l] and [2] This is achieved by mul¬ 
tiplying the expected fraction of sources in a region of 
the phase-space with the number of sources expected to 
occur during the search. This gives the actual number 
of sources detected. As expected, Type la SNe are the 
most numerous contaminant with ^3 — 12 sources ex¬ 
pected to exhibit i — z > 0 mag. Strategies 3 and 4 
produce the fewest number of Type la SNe detections 
3 — 5), while Strategy 2 produces the most 12). 
For the Pan-STARRS fast transients, we find that the 
number of sources exhibiting i — z > 0 mag is ^ 3 
across all four strategies. Lastly, we expect 1 type 
.la and miscellaneous fast transients to exhibit i — z > 0 
mag. If we additionally require that the rise time is < 4 
days the number of detected Type la SNe is reduced to 
~ 1 — 7. This does not affect the number of detected Pan- 
STARRS fast transients. Ultimately, we do not expect 
the lack of a pre-existing template to significantly hinder 
our ability to separate kilonovae from contaminants on 
the basis of their i — z color and timescale, if the source 


is observed during the rise to peak. 

Figure [l5| shows the rise time-raj slice of our phase- 
space. We first note the wide range of values for raj that 
appear in our simulations. This arises from considering 
contaminants over a wide redshift range or in the case 
of the kilonovae, over a range of model parameters. We 
further note that there does not appear to be an obvious 
cut utilizing raj that helps to distinguish kilonovae from 
the contaminant population. However, we can separate 
kilonovae from the contaminant populations by utilizing 
a cut on timescale as all of the kilonovae exhibit a rise 
time of < 4 days. Therefore, if we require that sources 
exhibit a rise time of < 4 days, then we expect to elim¬ 


inate ~ 100 Type la SNe, ~ 1 Pan-STARRS fast tran¬ 
sient, and < 1 type la and miscellaneous fast transients, 
without rejecting any kilonovae. 

Despite the lack of a clear separation between the kilo¬ 
novae and contaminant populations in raj, the require¬ 
ment that sources exhibit a raj greater than 5cr in the 
difference images does reduce the overall number of con¬ 
taminants detected. Utilizing this criterion, the detec¬ 
tion rates for Type la SNe and Pan-STARRS fast tran¬ 
sients (the two most numerous contaminants) on the rise 
is ~ 1% for both populations, down from ~ 60% and 
~ 25%, respectively, in the case of a pre-existing tem¬ 
plate. The downside of this approach is that ~ 10 — 15% 
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Fig. 16. — Same as Figure [14] but showing the decline time-color slice of the phase-space. The kilonovae are not always clearly separated 
from the contaminant population due to the imprecise color measurements. However, contaminants (with the exception of the miscellaneous 
fast transients) do stand out from kilonovae on the basis of their long decline times (> 6 days, the maximal length given our search duration). 


of kilonovae are also lost. These sources can be recovered 
with a late-time template, but they will not be detected 
in real time. 

In Figure [T6| we plot the decline time-color slice of our 
phase-space. We find that there is no longer a clear sep¬ 
aration in color between the kilonovae and contaminant 
populations. This is true for all contaminants and across 
all observing strategies. This effect is the result of the 
apparent slow evolution of the light curve during the de¬ 
cline phase relative to the template flux. This can occur 
if the observations straddle the light curve peak or if 
the source simply evolves slower after peak. However, 
we can separate kilonovae from the contaminant popula¬ 
tions by utilizing a cut on timescale as all of the contam¬ 
inants predominantly exhibit a decline time of > 6 days 
(the maximal value given the duration of our search). 
Therefore, if we require that sources exhibit a decline 
time of < 5 days, then we expect no Type la SNe, type 
.la and Pan-STARRS fast transients, and <C 1 miscel¬ 
laneous fast transient to be detected, but this cut will 
reject ^ 15 — 20% of kilonovae. 

Figure pT| shows the decline time-mj slice of the phase- 
space. As m the case of sources detected on the rise, 
there is no clear cut on that will help to separate the 
kilonovae from any contaminants. Furthermore, a cut on 
decline time (e.g., > 6 days) does not yield any new infor¬ 


mation beyond what is gained from the decline time-color 
slice of the phase-space. However, as before, utilizing 
a detection criterion on ra* does reduce the total num¬ 
ber of contaminants detected. This effect is enhanced 
for sources such as Type la SNe and the Pan-STARRS 
fast transients, as they show slower post-peak evolution 
in their light curves. Here, the detection rates for these 
sources, when seen on the decline, is <C 1 % in both cases. 
This is down from ~ 20% and ^ 15%, respectively, in the 
case of a pre-existing template. However, the detection 
criterion on m* reduces the kilonovae detection rates by 
-5%. 

We now summarize the efficacy of cuts based on i — z 
color, timescale, and change in flux for separating kilono¬ 
vae from the contaminants when a pre-existing template 
is not available. For sources observed during the rise to 
peak, requiring i — z > 0 mag eliminates ~ 98% of de¬ 
tected Type la SNe, ~ 91% of detected type .la, ~ 100% 
of detected miscellaneous fast transients, and ~ 56% of 
detected Pan-STARRS fast transients while eliminating 
< 0.1% of detected kilonovae. If we additionally require 
that the rise time is < 4 days, then we eliminate ^ 99% of 
detected Type la SNe, ~ 94% of detected .la SNe, and 
^ 59% of detected fast Pan-STARRS transients, with¬ 
out altering the fraction of kilonovae eliminated. In this 
scenario, an additional cut on ra* does not yield any ad- 
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Fig. 17. — Same as Figure [I5| but showing the decline time-ra* slice of th e phase-space. As in Figure [l5| the contaminants are not clearly 
separated from the kilonovae population in m*. However, as in Figure [l6| the long decline times can be used to distinguish kilonovae from 
the contaminant populations. 


ditional information. 

For sources observed only in decline, a cut of i — z > 0 
mag eliminates ^ 50% of detected Type la SNe, ^ 99% of 
detected type .la, and ^ 100 % of detected miscellaneous 
fast transients. The fraction of detected Pan-STARRS 
fast transients eliminated depends on the observing strat¬ 
egy; for Strategies 1 and 3 it is essentially zero, while for 
Strategies 2 and 4 it is ~ 11% and ~ 37%, respectively. 
The reason for this variation is the fact that the different 
observing cadences are better matched to specific subsets 
of the population. The cut of i — z > 0 mag color also 
eliminates ^ 15% of detected kilonovae. If we addition¬ 
ally require that the decline time is < 5 days then we 
remove ^ 100 % of the contaminant population but also 
eliminate ^ 60 — 80% of detected kilonovae. No addi¬ 
tional information is gained from a cut on m*. There¬ 
fore, cuts on i — z color and timescale are less effective 
when the source is observed on the decline only. This 
highlights the need for the rapid triggering of follow-up 
observations. 

Despite these challenges, the best approach in the ab¬ 
sence of a pre-ex istin g template is unchanged from the 
results of Section |4.2| The optimal follow-up cadence is 
still Strategy 1 , which allows good sampling of the light 
curve while maximizing the area that can be covered in 
a single night of observations. This strategy is also not 


subject to the dip in detection efficiency as a function 
of start time, seen in Strategies 2 and 4. However, we 
note that the previously recommended depths of 24 AB 
mag in z-band and 23 AB mag in z -band are only ex¬ 
pected to detect ^ 50 — 75% of kilonovae in this scenario 
(Figure 12 ). The z — z color still provides a strong con¬ 
straint for separating kilonovae from contaminants but it 
is most reliable when the kilonova can be observed dur¬ 
ing the rise to peak. If this is not feasible due to other 
observational constraints, then it is best to resort to ob¬ 
taining template images once the kilonova has faded in 
the optical (> 10 days). 

6 . EFFECTS OF ALTERNATIVE KILONOVAE MODELS 


So far, we have only considered kilonova models af¬ 
fected by the large opacities of r-process heavy nuclei. 
We now consider the effects of modifications to the stan¬ 
dard kilonova models and their impact on our detectabil¬ 
ity results. The key issue in determining the nature of 
kilonova emission is the electron fraction, Y e . The value 
of Y e in the dynamical ejecta and from accretion disk 
winds is low (Y e <0.2) which leads to the production 
of r-process nuclei and strong suppression of the optical 
emission. However, if the electron fraction of the ejecta is 
above the critical value of Y e ~ 0.25, r-process nucleosyn¬ 
thesis is not able to produce heavy elements (particularly 
the lanthanides, Kasen et al. 2014). The resulting opac- 
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ities will be closer to those from the Fe-peak elements, 
leading to a bluer transient. Here we investigate spec¬ 
ulative scenarios that could produce a “blue” kilonova 
component and assess their impact on r-band detectabil¬ 
ity. 

6.1. Neutrino Driven Winds 

Ejecta can be produced by the outflows from a post¬ 
merger accretion disk (Dessart et al. 2009, Fernandez & 
Metzger 2013, Fernandez et al. 2014). Fernandez & Met¬ 
zger (2013) performed numerical two-dimensional hydro- 
dynamical simulations of accretion disks and found that 
^10% of the initial disk mass is ejected over timescales of 
~ 1 s. These ejecta are a neutron-rich wind with an elec¬ 
tron fraction of Y e ~ 0.2, producing conditions that are 
favorable for r-process nucleosynthesis. This type of out¬ 
flow will therefore not drastically alter the observational 
properties of the kilonova. 

Recent work by Fernandez et al. (2014) found that 
the amount of ejected mass is a function of the spin 
of the remnant black hole, with more rapidly spinning 
black holes leading to higher ejecta masses by a factor 
of several. Furthermore, they found that Y e increased 
monotonically with black hole spin due to an increased 
early-time neutrino luminosity from a rapidly spinning 
black hole. They found that in the most extreme cases 
(a > 0.8), the system could eject ~ 10 -5 — 3 x 10 -3 M 0 
of material with Y e > 0.25. Specifically, an ejecta mass of 
^ 10 -3 M 0 will produce a blue component in the kilono- 
vae light curve with uL v ~ 7.5 x 10 4 ° erg s -1 at ^ 6500 
A(r ~ 22.5 mag at 200 Mpc) and a peak time of < 1 
day (Kasen et al. 2014). This is significantly brighter in 
r-band than the kilonovae models discussed previously, 
but the short timescales make observing this component 
challenging. 

6.2. A Long-Lived Hyper-massive Neutron Star 

Another modification can be caused by a surviving 
hyper-massive neutron star (HMNS). The HMNS is a 
potential result of NS-NS mergers, but in simulations it 
often collapses rapidly (t < 100 ms) to form a black hole 
(Sekiguchi et al. 2011). The prompt formation of a black 
hole is an implicit assumption in the kilonova models dis¬ 
cussed so far in this paper. However, if the total mass 
of the binary is less than a critical value, the HMNS 
can survive for a longer period of time (Bauswein et al. 
2013a, Kaplan et al. 2013). In this scenario, the resulting 
neutrino luminosity will raise the electron fraction of the 
disk wind ejecta and will suppress r-process nucleosyn¬ 
thesis. Thus, the lanthanide-free disk wind material will 
produce a brighter, bluer, and shorter lived component in 
the kilonova light curve (e.g., Metzger & Piro 2014, Met¬ 
zger & Fernandez 2014, and references therein). This 
emissions component will have a strong viewing angle 
dependence, with the polar regions being more strongly 
irradiated by neutrinos than the equatorial regions which 
will remain obscured by lanthanide-rich material. Con¬ 
sequently, the blue early component is more pronounced 
when the merger is observed along the polar axis. 

Metzger & Fernandez (2014) found that the luminosity 
and timescale of this initial blue component are strongly 
dependent on the lifetime of the HMNS, with typical 
values (^ 100 ms) leading to peak luminosity that is a 


factor of a few times brighter than the red kilonova com¬ 
ponent, with a timescale of a few days for emission along 
the polar direction. The light curves of this early-time 
component were computed in detail using wavelength- 
dependent radiative transfer calculations by Kasen et al. 
(2014). They found that for an HMNS with a lifetime of 
^100 ms, the initial blue component will have a peak 
luminosity of vL u ^ 1.3 x 10 41 ergs s -1 at ~ 6500 A 
(r ~ 22 mag at 200 Mpc) and a peak time of ^ 2 days. 
Kasen et al. (2014) also compared the blue component 
to the r-process powered emission from a 10 -2 M 0 torus 
of dynamical ejecta. They found that the blue compo¬ 
nent will be visible for polar viewing angles (0 < 15 deg). 
The fraction of sources expected to be observed at this 
viewing angle is ^ 7 0 2 ~ 0.48 (Metzger & Berger 2012). 
However, at these viewing angles the kilonova emission 
will be overwhelmed by the optical afterglow of an on- 
axis SGRB that may accompany the merger. Further¬ 
more, if the source is observed at large viewing angles 
(6 > 15 deg) then the emission will become obscured by 
the dynamical torus and the r-band emission will then 
peak at > 23.7 mag (~ 2 mag fainter than the polar case 
for a torus mass of 10“ 2 M 0 ). 



Fig. 18. — r-band ligh t cur ves from the alternative kilonova mod¬ 
els discussed in Sections |6.2| and |6.3| The r-band light curve for the 
neutron-rich dynamical ejecta with M e j = 10 -2 Mq and /3 e j = 0.2 
is plotted for comparison. The alternative models produce light 
curves which can be significantly brighter and bluer than the dy¬ 
namical ejecta but peak at earlier times. 


6.3. A Neutron Precursor 

Lastly, we consider the case in which a small frac¬ 
tion of the neutron-rich ejecta is expanding so rapidly 
that it avoids the initial phase of r-process nucleosyn¬ 
thesis. This effect has been seen in several general rel¬ 
ativistic simulations (Bauswein et al. 2013b, Goriely et 
al. 2014, Just et al. 2014). In particular, Bauswein et 
al. (2013b) found that for a small fraction of the ejecta 
(^ 10 -5 — 10 -4 Mq), the expansion timescale was suffi¬ 
ciently rapid (< 5 ms) that they evade neutron capture. 
These neutrons will undergo /3-decay and release energy 
into the ejecta. Given that the photon diffusion time is 
on the order of a few hours in the outer layers of the 
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Fig. 19. — Contour plots for t he r -band detectability of alternative kilonovae models. The top row shows the contours for the dynamical 
ejecta (reproduced from Figure [4T] >. The middle row shows contours for a long-lived HMNS with a lifetime of ~ 100 ms, surrounded by a 
dynamical tours with a mass of 10 -2 Mq of neutron-rich material observed at a viewing angle of 9 ~ 18 deg as modeled by Kasen et al. 
(2014). The bottom row shows the effect of an early-time neutron precursor as modeled by Metzger et al. (2014). 


ejecta, a significant fraction of this energy will manage 
to produce radiation. The light curves of this “neutron 
precursor” were modeled by Metzger et al. (2014). They 
found that the /3-decay from the population of free neu¬ 
trons at the leading edge of 10 -4 M 0 of ejecta can pro¬ 
duce a transient with a timescale of ~ 1 — 2 hours and a 
peak r-band magnitude of ~ 22.2 mag at 200 Mpc. 


6.4. Implications for r-band detectability 


We now investigate the effect of these speculative mod¬ 
els on the early-time detectability of kilonovae using r- 
band observations. We construct the light curves for the 
long-lived HMNS from the time-resolved spectra mod¬ 
eled by Kasen et a l. (2014) using the procedure out¬ 
lined in Section 2.1 These spectra are for a HMNS life¬ 
time of 100 ms with a dynamical torus with 10 -2 M 0 
of neutron-rich material observed at a viewing angle of 
cos# = 0.95 (0 ~ 18 deg). For the neutron precur¬ 
sor, we use the r-band light curve provided by Metzger 
et al. (2014), for 10 -4 M 0 of free neutrons with an 
electron fraction of Y e « 0.05 and opacity of k « 30 
cm 2 g -1 . These light curves are plotted in Figure fl8| 
along with those from the dynamical ejecta-only model 
for M e j = 10 -2 M 0 and /3 e j = 0.2 used in Sections [il and 

El 


We r epea t the detectability simulations discussed in 
Section |4.1| for the two speculative light curves in Fig¬ 
ure [18] which are co-added in flux space to the light curve 
for the dynamical ejecta from Barnes & Kasen (2013). 
The general simulation procedure and criteria for d etec- 
tion are identical to those described in Section [4711 The 
contours are plotted in Figure [l9j The r-band contours 
for the dynamical ejecta only (Figure [4| are reproduced 
for comparison. 


Inspecting the case of a long-lived HMNS, we first note 
that the typical depth required to detect 50% (95%) of 
events is 23.4 (24) mag across all four observing scenar¬ 
ios. This is shallower than the depth required for the dy¬ 
namical ejecta alone, but still within reach of the larger 
telescopes. When a slower cadence strategy is utilized 
(e.g., Strategy 2 which involves observations every other 
night after the first observation) the 95% detection cri¬ 
terion drops to 24.7 mag for observations starting more 
than six hours after the GW trigger. This is due to the 
rapid decline of the light curve after peak. Observations 
must then be deeper at late times to obtain the necessary 
three epochs for a detection. If a rapid cadence is utilized 
(e.g., Strategies 3 and 4), then the 50% detection fraction 
contour is at 23 mag for observations within six hours of 
the GW trigger. The other contours are unaffected. 

For the neutron precursor model, the typical depth re¬ 
quired to detect 50% (95%) of events is 24.7 (25.9) mag 
for observing Strategies 1 and 2, comparable to those re¬ 
quired for the dynamical ejecta alone because the precur¬ 
sor is short-lived. When Strategies 3 and 4 are utilized, 
the 50% detection fraction contour begins at 23.6 mag 
for observations starting 3 hours after the GW trigger. 
It then monotonically declines to 24.7 mag by 24 hours. 
This is again because a neutron precursor is short-lived, 
and only early time observations benefit from its high 
luminosity. 

Contour plots showing the detecti on o f a rise to peak 
are shown in Figure |20| As in Section [4T| a rise is defined 
as any observed brightening between epochs. We note 
that for Strategies 1 and 2, both of which involve a single 
observation on the first night, the source is always seen 
in decline for all three models. In the case of a long-lived 
HMNS brightening is observed in 50% (95%) of events for 
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Fig. 20.— Same as Figure [I9| but showing the fraction of events for which a brightening is detected. Compared to the case of dynamical 
ejecta only, it is significantly more challenging to detect brightening in these alternative kilonova models due to their earlier peak times 
and rapid evolution. 


observations starting < 5 (< 3) hours of the GW trigger 
and reaching a depth of ~ 23 (24) mag. For a neutron 
precursor, brightening in the source is never observed. 
This is because the source peaks in ^ 2 hours, while a 
reasonable start time is longer. 

Ultimately, the case of a long-lived HMNS irradiating 
the ejecta shows the strongest potential impact on kilo- 
nova r-band detectability. However, this component has 
a strong dependence on viewing angle and thus will not 
always be observed. In addition, this component is still 
faint, requiring observations that achieve a depth of « 24 
mag in r-band to detect 95% of the sources. This scenario 
is highly speculative and requires further study of the 
neutron star equation of state and post-merger behavior 
to assess the likelihood that the post-merger remnant will 
avoid the prompt collapse to a black hole. The scenario 
involving a neutron precursor does not drastically impact 
r-band detectability unless an observing strategies with 
rapid intra-night cadence is employed (i.e., a baseline of 
~ 1 hr). 

7. FOLLOW-UP STRATEGIES 

Motivated by the results of our simulated observations, 
we now assess the most effective methodologies for con¬ 
ducting efficient GW follow-up observations. There are 
two primary issues to address: The detection of a kilo- 
novae, and its identification against a population of con¬ 
taminant sources. Our simulations show that, in the case 
of a pre-existing template, observations which achieve a 
depth of i « 24 mag and z « 23 mag lead to the detection 
of about 95% of kilonovae. Furthermore, observations 
should commence within 12 hours of the GW trigger to 
maximize the likelihood of detecting the source during 
the rise to peak. We also show that cuts made on i — z 
color, the rise time (if available), and A m z are strong 


discriminants for separating a kilonova from the contam¬ 
inant population. Quantitatively, with i — z > 0 mag, 
a rise time of < 4 days, and \Am z \ > 0.4 mag, we can 
eliminate a significant fraction (> 95%) of contaminants 
without removing any kilonovae. If the source is seen in 
decline then requiring i—z > 0 mag and | A m z | > 0.4 mag 
are still strong cuts, eliminating ~ 60% of contaminants. 
If we also require a decline time of < 5 days then we 
can eliminate > 90% of contaminants, however this cut 
eliminates nearly half of the kilonova population. These 
results do not depend strongly on the choice of cadence 
(i.e., there is no benefit from intra-night observations). 
Therefore, follow-up observations should include a single 
visit per filter per night for at least the first few nights 
following a trigger. 

Taking these constraints into account; when developing 
the choice of observation strategy there are three ques¬ 
tions that must be considered: 

1. How large of an area must be covered by the ob¬ 
servations? 

2. How long is the region of interest observable? 

3. How soon after the GW trigger can observations 
commence? 

The first and second points are the most important 
considerations as the sizable localization regions from the 
early science-runs of Advanced LIGO present the largest 
observational challenge and have the most pronounced 
impact on the quality of data that can be obtained. The 
key issue here is whether or not there will be sufficient 
time to observe the entire GW localization region on a 
nightly basis. A successful GW follow-up campaign must 
be able to observe a significant fraction of the GW lo¬ 
calization region in a single night to obtain the optimal 
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TABLE 4 

Follow-up Capabilities of Wide-Field Telescopes 


Instrument 

FoV 

(deg 2 ) 

Etendue 
(m 2 deg 2 ) 

r- 

Exp. 

(s) 

■band 

Rate 

(deg 2 hr -1 ) 

Exp. 

(s) 

z-band 

Rate 

(deg 2 hr -1 ) 

Exp. 

(<0 

z-band 

Rate 

(deg 2 hr -1 ) 

z + z Rate 
(deg 2 hr -1 ) 

DECam a 

3 

37.7 

4200 

2.5 

250 

43 

150 

72 

27 

ZTF b 

47 

53.2 

9 x 10 5 

0.2 

1 .6xl0 4 

10 

— 


— 

Subaru HSC C 

1.8 

90.5 

520 

12 

55 

115 

45 

1 15 

65 

Pan-STARRS2 d 

14 

35.6 

1.4 x 10 4 

4 

700 

72 

350 

144 

30 

LSST e 

9.6 

315 

390 

88 

28 

1235 

19 

1820 

735 


Note. — Expected rates of sky coverage in deg 2 hr -1 for several existing and future wide-field facilities at a depth of r « 26 mag, 
i ~ 24 mag and z ~ 23 mag (i.e., the depths necessary to achieve a 95% detection rate for kilonovae). We note that for HSC and LSST we 
include overhead times of 20s and 2s, respectively, as they are comparable to the exposure times. References: 


a DECam ETC, http://www.ctio.noao.edu/noao/content/Exposure-Time-Calculator-ETC-0 
b Nissanke et al.(2013), Bellm (2014), Kasliwal & Nissanke (2014) 
c HSC ETC, http://www.naoj.org/cgi-bin/img_etc.cgi 
d Pan-STARRS DRM, http://pan-starrs.ifa.hawaii.edu/ 

e LSST ETC, http://dls.physics.ucdavis.edu:8080/etc4_3work/servlets/LsstEtc.html 


cadence. The third question impacts the likelihood that 
we will be able to observe a kilonovae bef ore p eak. This 
is important as we have shown in Sections |4. 2 1 and |5| that 
it is easier to isolate a kilonova from the contaminant 
population if it is observed during the rise. Additionally, 
observing a kilonova at early times is important both for 
constraining the ejecta parameters and for testing the 
alternative kilonovae models outlined in Section [6j Ulti¬ 
mately, our ability to address the three questions above 
depends heavily on the facilities in use. 

To investigate the impact of the choice of observing 
facility on GW follow-up we determine the expected per¬ 
formance of several current and future wide-held tele¬ 
scopes. Specifically, we consider the performance of DE¬ 
Cam, ZTF, Subaru Hyper Suprime-Cam (HSC), Pan- 
STARRS2, and LSST. We take into account the field-of- 
view (FoV) per pointing, the required integration time to 
achieve a depth of r ~ 26, i ~ 24, and z ~ 23 mag, and 
consequently how much area can be covered per hour of 
observing time. We also compute the sky coverage rate 
when observing in i- and z -band to obtain color data. 
These values are given in Table |4j 

We first look to currently available facilities (e.g., DE¬ 
Cam, HSC, and Pan-STARRS2). In the Northern hemi¬ 
sphere the two available telescopes are Pan-STARRS2 
and HSC. Pan-STARRS2 offers an impressive FoV 14 
deg 2 ) allowing a ^ 100 deg 2 region to be imaged in just 8 
pointings. However, HSC is the clear front-runner as the 
immense collecting area allows observations to achieve 
the target depths in much shorter integration times than 
Pan-STARRS2. Thus, despite the smaller FoV of HSC, 
the sky coverage rate for i- and 2 -band data is faster than 
Pan-STARRS2 by over a factor of two. If we consider a 
typical GW localization region 100 deg 2 ) then Pan- 
STARRS2 can image this region to the required i- and 
z -band depths in ^ 3.5 hours while achieving this goal 
with HSC will require only ^1.5 hours. In the Southern 
hemisphere DECam is the ideal choice. The wide FoV 
and red-sensitive CCD of DECam allow it to quickly im¬ 
age a typical GW localization region. If we again consider 
a ~ 100 deg 2 localization region, then DECam can ob¬ 
tain i- and 2 -band data at the required depths in ^ 3.5 


hours. 

We now consider future wide-held telescopes. The next 
major instrument to become available will be ZTF. The 
primary advantage of ZTF is the extremely wide FoV 
47 deg 2 ), allowing even large localization regions to 
be imaged in just a handful of pointings. However, the 
longer integration times required to achieve the target 
depths means that the sky coverage rate of ZTF will 
be slower than DECam, Pan-STARRS2, and HSC in r- 
and Aband. Furthermore, the lack of z-band capabilities 
makes obtaining i — z color data impossible (observations 
in Tband alone may still be useful, see below). Looking 
further ahead, GW follow-up with LSST will be partic¬ 
ularly effective. The large collecting area and high sen¬ 
sitivity allow LSST to image a GW localization region 
with unparalleled rapidity. If we again consider a typical 
localization region (^ 100 deg 2 ), then ZTF will require 
~ 500 hours in r-band and and ~ 10 hours in i-band to 
image the entire region to the required depths. By com¬ 
parison, LSST can image the entire localization region to 
the required depths in i- and z-band in just ~ 10 min¬ 
utes, highlighting the immense observing power of future 
facilities. However, LSST is not expected to begin sci¬ 
ence operations until ~ 2022 when GW localizations will 
be much improved. 

While GW follow-up observations should strive to ob¬ 
tain i — z color data as outlined above, this may not be 
possible because either there is insufficient time to ob¬ 
serve the localization region in i + z or because z-band 
observations are not feasible (e.g., ZTF). In this scenario, 
data may be taken in a single filter, the identification of 
a kilonova becomes difficult. We dedicate the rest of this 
section to a brief discussion of eliminating contaminants 
when only single filter observations are possible. 

The most populous contaminant to contend with is 
the Type la SNe. In the context of observations uti¬ 
lizing a single filter, these contaminants are most easily 
eliminated by utilizing their much slower rise and de¬ 
cline timescales. Therefore, the best method to reject 
these sources is to perform image subtraction using ob¬ 
servations that are separated by (e.g., > 1 day between 
epochs). Type la SNe are unlikely to be detected in these 
subtractions as they are unable to produce a statistically 
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significant change in flux between the two epochs. Qual¬ 
itatively, if we require a source to exhibit a change in 
flux of at least \Am z \ > 0.4 mag and a rise time of < 4 
days, then we can eliminate 99% of detected Type la SNe. 
This corresponds to ~ 6 detected events remaining after 
the cuts have been applied. However, we also eliminate 
^ 15 — 35% of detected kilonova, with the higher frac¬ 
tion of sources eliminated in Strategies 3 and 4 where the 
rapid cadence of observations on the first night produces 
smaller values of \Am z \. For sources seen on the decline, 
if we require that the decline time is < 5 days and the 
change in flux is \Am z \ > 0.4, then we eliminate all of 
the detected Type la SNe, while rejecting only ~ 2% of 
detected kilonovae. In both scenarios (e.g., sources seen 
to rise or only in decline) rejected kilonovae can be re¬ 
covered using a late-time template image, but real time 
detection may be difficult. Still, despite their high ex¬ 
pected rates, Type la SNe will likely prove to be one of 
the easiest contaminants to eliminate in GW follow-up 
searches. 

We can also apply the above argument to Pan- 
STARRS fast transients that are seen only in decline. 
Qualitatively, if we require that sources exhibit a decline 
time of < 5 days and a change in flux of \Am z \ > 0.4, 
then we can eliminate ^ 80% of the detected Pan- 
STARRS fast transients, again while rejected ~ 2% of 
the detected kilonovae. While the cuts on decline time 
and | A m z | are not as effective here as in the case of Type 
la SNe, the lower rates 7 total expected during a typi¬ 
cal GW follow-up search) means that we expect only ^ 2 
detected events to survive. These remaining sources can 
then be followed-up and eliminated by utilizing pointed 
observations as described below. 

While inter-night variability is an effective tool for 
eliminating slowly evolving contaminants, it is not fea¬ 
sible for the rapidly-evolving contaminants (e.g., type 
.la, miscellaneous fast transients, and Pan-STARRS fast 
transients observed on the rise) as they exhibit timescales 
and changes in brightness comparable to those of kilono¬ 
vae. We can rule out type .la and miscellaneous fast 
transients as likely cont amina nts by utilizing the rate in¬ 
formation from Section [2. 101 We showed that the rates 
for such events are low, resulting in < 1 source expected 


on the timescale and in the area of a typical GW follow¬ 
up campaign. We further showed in Section |5j that if the 
rate information is convolved with the detection rates 
in our simulations, then we expect < 1 of these con¬ 
taminants to appear during the course of follow-up ob¬ 
servations. In this case, even considering the fact that 
the rates are poorly constrained, it is unlikely that these 
sources will lead to significant contamination. However, 
this rate argument does not apply to the Pan-STARRS 
fast transients as the rates of these events are sufficiently 
high 7 expected during follow-up of a typical GW er¬ 
ror region) that they are likely to represent significant 
contamination. In this case, the best option is to in¬ 
clude any such sources in a list of kilonova candidates 
for follow up with pointed observations as described be¬ 
low or wait until the source begins to decline and apply 
the inter-night variability arguments outlined above. 

La stly , we consider flaring stars (e.g., M-Dwarfs, Sec¬ 


tion 2.8) as such events are likely to become contami¬ 
nants m observing strategies that rely on a single epoch 
to identify sources. This is particularly true for observa¬ 


tions carried out in bluer bands (e.g., g or r) where the 
flares amplitudes are larger and the kilonova timescales 
are much shorter. Given that flares generally exhibit 
short timescales (minutes to hours, Berger et al. 2013) 
then if we require that a source be detected in at least 
two epochs separated by > 1 day we can efficiently elim¬ 
inate all flare contaminants. The primary concern here 
is that rapidly evolving kilonovae (e.g., in r band around 
peak) will also be eliminated, particularly if observations 
are unable to achieve a depth of r « 26 mag resulting in 
the kilonovae only being detected in one or two epochs. 
Additionally, in r-band, the quiescent counterpart may 
be difficult to identify in the template images. These 
complications can be avoided by carrying out follow-up 
observations in i- or z-band. 

Taking the above considerations into account, we con¬ 
clude that single filter observations, particularly those 
conducted in z-band, are still viable for kilonova identi¬ 
fication. We favor data taken in z-band as the shallower 
depths required to achieve a 95% detection rate (when 
compared to the depths required to achieve a compa¬ 
rable detection rate in r- or i-band) allow substantially 
more sky area to be covered in a night of observations. 
Furthermore, we have shown above that observations in 
z-band still allow robust contaminant rejection (although 
still less robust than when color data is available) with¬ 
out being susceptible to additional contamination from 
stellar flares as may be the case for observations in r- 
and i-band. We note that observational strategies that 
utilize a single filter are not likely to produce a single, 
obvious kilonova. This may be particularly true when 
a pre-existing template is unavailable, reducing the effi¬ 
cacy of our cuts on the timescale and change in bright¬ 
ness. However, the strategies outlined in this section 
can still be utilized to generate a list of kilonova candi¬ 
dates. In this scenario, once a list of < 10 candidates is 
generated, those sources can be further observed in the 
NIR and spectroscopically. NIR follow-up is crucial as 
kilonova are substantially redder than the other contam¬ 
inants considered, with the kilonova spectrum peaking at 
~ 1.5/im. Furthermore, spectroscopic observations will 
allow us to look for features consistent with the produc¬ 
tion of r-process elements as another means of kilonova 
identification. These observations will vastly improve the 
robustness of kilonova identification and will allow us to 
maximize the science gained from a GW event. 

The worst-case scenario is a GW localization region 
that is simply too large to observe at the required depth 
in the available time, even in a single filter. Here, there 
are two possible approaches. The first is to simply obtain 
shallower observations at a depth that allows the entire 
GW localization region to be imaged. This will be useful 
if the source is located closer than the fiducial distance of 
200 Mpc considered here, or if the ejecta mass is larger 
than the fiducial value of 0.01 M 0 . However, this is 
risky as shallower observations are unlikely to detect the 
typical event. Furthermore, even if a kilonova is detected, 
it may only appear in one or two epochs, making it more 
difficult to distinguish from stellar flares. A less risky 
approach is to carry out the observations at the required 
depth, but cover the highest probability portion of the 
localization region, thereby maximizing the likelihood of 
detecting the kilonova. 
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8. CONCLUSIONS 

The era of gravitational wave astronomy is fast ap¬ 
proaching, with ground-based GW interferometers ex¬ 
pected to achieve direct detections of the inspiral and 
merger of compact object binaries sometime during the 
next few years. If we wish the maximize the science gains 
of these detections then the rapid identification of an EM 
counterpart is essential. Observational and theoretical 
evidence suggests a host of potential EM counterparts 
spanning a wide range of timescales, luminosities, and 
wavelengths. Among these, the isotropic optical emis¬ 
sion powered by the radioactive decay of r-process el¬ 
ements synthesized in the merger ejecta has proven to 
be a particularly enticing target. A major challenge for 
follow-up observations is the large sky area that needs 
to be observed 100 deg 2 ) combined with a substan¬ 
tial depth to detect a typical event. Consequently, op¬ 
tical searches for kilonovae will be plagued by contami¬ 
nation from other common transient events (e.g., Type 
la SNe) as well as from other rapidly evolving transients 
(e.g., see Section [2|. We have performed comprehen¬ 
sive simulations designed to probe the range of depths 
and cadences that may be involved in a GW follow-up 
campaign to determine the set of observations that will 
allow kilonove detections coupled with efficient rejection 
of contaminants. We find that: 


• In the case where a pre-existing templ ate i mage 
of the target region is available (Section 4.1), the 
optimal strategy for kilonova detection is to im¬ 
age the region to i « 24 and z & 23 mag with a 
single visit per night. We have shown that these 
limits are necessary to achieve a 95% kilonova de¬ 
tection rate. Furthermore, observations should ide¬ 
ally commence within 12 hours of the GW trigger 
to facilitate the detection of a rise to peak in the 
light curve. 


• For sources observed on the rise, cuts on i — z color 
and rise time are effective at distinguishing kilono¬ 
vae from the contaminant population (Section |4.2| ). 
Specifically, with i — z > 0 mag and a rise time of 
< 4 days, it is possible to eliminate > 95% of con¬ 
taminants without rejecting any kilonovae. Apply¬ 
ing the same cuts to sources seen only in decline is 
less effective. If we require i — z > 0 mag and a 
decline time of < 5 days then we again eliminate 
> 95% of the contaminant population, but lose 
nearly half of the kilonovae. These rejected kilo¬ 
novae can be recovered using a late-time template 
(> 10 days), but real-time kilonova identification 
will be hampered. 


• In the case where a pre-existing template is not 
available (Section^, we find that there is no effect 
on the optimal caaence, but observations must go 
^ 1 magnitude deeper than in the pre-existing tem¬ 
plate case to achieve comparable detection rates. 
The detection of a rise is also more difficult since 
the first observation must be used as a template 
image, reducing the amount of early-time data. 

• In the case where a pre-existing template template 
is not available, cuts on i — z color and timescale are 


still effective for separating kilonovae from the con¬ 
taminant population for sources observed during 
the rise to peak (Section [5|. However, for sources 
seen only in decline, these cuts are significantly less 
effective. We do note that by utilizing a detection 
criterion based on the change in brightness in an 
observed source we can reduce the number of con¬ 
taminant sources detected. However, this also re¬ 
duces the kilonovae detection rates. 

• For the alternative kilonova models (Section we 
find that r-band detectability in the speculative 
case of a long-lives HMNS producing an early-time 
blue component may improve, but the short-lived 
nature of this component combined with the strong 
dependance on viewing angle means that it is not 
likely to be observed. If the kilono va is accompa¬ 
nied by a neutron precursor (Sect ion [6~3] ) , then this 
component will appear brighter than a kilonova 
(r « 22.2 mag), but the extremely short timescale 
(< 2 hrs) makes it unlikely to appear in follow-up 
observations. 


• We examine the potential of current and future 
wide-held telescope facilities to achieve the observa¬ 
tional goals outlined above (Section^. For current 
facilities, we find that Subaru HSC is the best op¬ 
tion for localization regions in the Northern hemi¬ 
sphere while DECam is best suited to GW follow¬ 
up in the Southern hemisphere. Both telescopes 
can easily image an entire GW localization region 
to the required depths in i- and z-band at the re¬ 
quired cadence. In the future, GW follow-up with 
LSST will be particularly effective and rapid. 


• If data in both i- and z-band cannot be obtained 
due to observational constraints, then an attempt 
at kilonova identification can be made using data 
from a single filter, ideally z-band (Section [7 ). We 
found that without color data it is still possible to 
eliminate the most numerous contaminants (e.g., 
Type la SNe and Pan-STARRS fast transients) by 
utilizing their slow light curve evolution when com¬ 
pared to kilonovae. The rapidly evolving contam¬ 
inants (e.g., type .la and miscellaneous fast tran¬ 
sients) are not likely to be a significant contaminant 
on the basis of their low expected rates. Lastly, M- 
dwarf stellar flares are easily rejected by requiring 
that the source appear in at least two epochs or by 
identifying a quiescent counterpart in the template 
image. 


The next step of this analysis is to test the method¬ 
ologies and conclusions of this work against a set of real 
data. As discussed in Section [l] we have obtained DE¬ 
Cam i- and z -band data with the cadence and depth op¬ 
timized for GW follow-up over an area of rsj 70 deg 2 . We 
will confront the results of this paper by using subsets of 
these data to mimic the four observing strategies as well 
as the effects of having or missing a pre-existing template 
(utilizing follow-up observations that we obtained about 
two weeks after the main survey). This will facilitate a 
better observational understanding of the true contami¬ 
nant population on timescales relevant to GW follow-up 
and further refine any proposed strategies. 
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